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1 Introduction 



The molecular structure of biological and biomimetic membranes is provided by 
bilayers of amphiphilic molecules such as lipids and proteins. Lipid bilayers are 
rather thin with a thickness between 4 and 5 nm. At physiological temperatures, 
these membranes are fluid in the sense that they have no in-plane shear modulus. 
On the nanometer scale, the molecules have no fixed neighborhoods since two 
adjacent molecules can easily exchange their position within the membrane. 
This local exchange leads to rapid lateral diffusion of the molecules along the 
membrane. In this way, a typical lipid molecule is displaced by about one 
micrometer in one second. Therefore, biomembranesj^represent multicomponent 
liquids in two dimensions. 

1.1 Membrane adhesion via sticker molecules 

Now, consider such a multicomponent membrane that forms a giant vesicle 
which adheres to another membrane or, more generally, to a substrate surface. 
Within the contact area of the vesicle, the membrane molecules can still diffuse, 
at least to some extent, which implies that the different molecular components 
may attain many different spatial patterns. In particular, some of these com- 
ponents may form clusters or extended domains within the contact area. In 
this article, we will review recent work on these patterns and their formation 
processes. 

We will focus on situations in which the membrane contains certain types 
of membrane-anchored molecules, so-called stickers and repellers. [1-6] Sticker 
molecules mediate attractive interactions between the membrane and the adja- 
cent surface whereas repeller molecules act as repulsive spacers between these 
two surfaces. Within the contact area, sticker molecules have a strong tendency 
to form clusters or intramembrane domains. In fact, there are several distinct 
mechanism for these pattern formation processes as explained in the main part 
of this review. 

There are many different types of stickers which can vary greatly in their 
size. The smallest stickers are presumably charged head groups of lipids which 
are attracted towards an oppositely charged surface. Likewise, lipids may have 
large sticky head groups containing, e.g., polysaccharides. Much larger stick- 
ers mediate the specific adhesion of biological membranes which governs both 
cell-cell adhesion and the signalling between cells of the immune system. [TJ [5] 
Cell adhesion molecules, which govern the binding of two cells, and receptor 
molecules, which are expressed on the surface of immune cells, are usually rel- 
atively stiff rodlike molecules which are anchored in the bilayer membrane and 
have a linear extension of 10 - 30 nm. 

^Here and below, the term 'biomembrane' is used as an abbreviation for 'biological or 
biomimetic membrane'. It is always assumed implicitly that biomembranes have two charac- 
teristic features: (i) they contain several molecular components and (ii) they are in a fluid 
phase. 
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Presumably the first theoretical models for the adhesion of membranes via 
mobile stickers were introduced by Bell et al. in Refs. [TU] as reviewed in 
[TT] and by Evans in Rcfs. [T2j[13]. In these models, the membrane is divided 
up into a bound segment, which represents the contact area, and an unbound 
segment which acts as a reservoir for the mobile sticker molecules which can 
diffuse in and out of the contact area. The repellers are taken to be more or 
less immobile and the separation of the two interacting surfaces is taken to be 
constant within the contact area. If the repeller molecules are longer than the 
sticker molecules, the competing action of these two types of molecules should 
lead to a modulation of this separation: sticker-rich membrane segments should 
have a relatively small separation whereas repeller-rich segments should have a 
relatively large one as discussed in the context of gap junctions O [T^ [T5] . 

Since the membrane-anchored stickers gain energy when they enter the con- 
tact area and bind to the second surface, these molecules will be 'recruited' 
by this surface and, thus, will be enriched within the contact area. Therefore, 
the bound and unbound segments of multicomponent membranes will usually 
differ in their molecular compositions [17||^ In the following, we will focus on 
the bound membrane segment, i.e., on the contact area of the multicomponent 
membrane, and view its unbound segment as a reservoir for the sticker and 
repeller molecules. It is then convenient to use a grand-canonical description 
and to replace this reservoir by chemical potentials for the sticker and repeller 
molecules. In this way, one arrives at lattice gas models on flexible surfaces. 
[U 0] These lattice gas models provide a general theoretical framework by 
which one can study the interplay of membrane adhesion and domain formation 
in a systematic manner. 

The problem of adhesion-induced domain formation has also been addressed 
using somewhat different theoretical models in Refs. [T^linillDJUT]. In addition, 
this process has also been observed experimentally both in biomimetic and in 
biological membranes. The next two subsections contain a short summary of 
these observations. 

1.2 Stickers and repellers in biomimetic membranes 

Adhesion-induced lateral phase separation into domains with small and large 
membrane separations has been found to occur in several biomimetic systems. 
The formation of blisters has been observed in membranes containing cationic 
lipids in contact with a negatively charged surface [22 , and between membranes 
containing both negatively and positively charged lipids 23J. The coexistence 
of tightly and weakly bound membrane segments has been found for mem- 
branes with biotinylated lipids bound to another biotinylatcd surface via strep- 
tavidin [24 , membranes with homophilic csA-receptors from the slime mold 
Dictyostelium discoideum [25j . and membranes containing specific ligands of in- 
tegrin molecules adsorbed on a substrate [26]. Attractive membrane- mediated 

^ This difference in composition has apparently been overlooked by Bell since he says on 
page 242 of Ref. [11) "that at equilibrium the number of receptors per unit area will be the 
same in the contact area and outside of it." 
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interactions between bound csA-receptors of adhering vesicles have been inferred 
from membrane tension jumps induced by the micropipet aspiration technique 
[27] . In addition to the receptors, the membranes studied in [211 [25l [261 127| also 

contain repulsive lipopolymers to prevent non-specific adhesion. 

1.3 Specific adhesion of biological membranes 

The adhesion of cells plays a key role in important biological processes such 
as tissue development and immune response. The highly selective interactions 
leading to cell adhesion are mediated by a variety of specific receptors which 
are embedded in the cell membranes. Two prominent examples for domain for- 
mation within the contact area of biomembranes are provided by focal contacts 
[71 I28j . which are formed during cell adhesion at substrate surfaces, and by 
the so-called 'immunological synapses' which are formed during the adhesion 
of helper T cells and antigen-presenting cells as part of the immune response. 
Within the contact area of these two cells, several groups have recently observed 
distinct domain patterns of membrane-anchored receptors and ligands. 

The antigen-presenting cells (APCs) display foreign peptide fragments on 
their surfaces. These peptide fragments are presented by MHC molecules on 
the APC surfaces, and recognized by the highly specific T cell receptors (TCR). 
At the cell-cell contact zone, the bound receptor- ligand pairs are arranged in 
characteristic supramolccular patterns [29l [301 ISl jSH [33] , for reviews see [SH 
|331|3Sll371|3Hll2ni • The final, 'mature' pattern of an adhering T cell consists of 
a central domain in which the TCRs are bound to the MHC-peptides (MHCp), 
surrounded by a ring-shaped domain in which the integrin receptors LFA-1 
of the T cell are bound to their ligands ICAM-1 of the APC. Intriguingly, 
the characteristic intermediate pattern formed earlier during T cell adhesion is 
inverted, with a TCR/MHCp ring surrounding a central LFA-1 /ICAM-1 domain 
in the contact zone [3011101133]. This pattern inversion has been first observed for 
T cells adhering to a supported lipid bilayer with embedded MHCp and ICAM-1, 
[3ni[10] , more recently also in a cell-cell system [33] . A significantly different type 
of pattern evolution has been found for immature T cells or thymozytes, which 
form multifocal synapses with several nearly circular clusters of TCR/MHC- 
peptide complexes in the contact zone [411 142] . 

1.4 Outline of review 

This review article is organized as follows. In section [2] we describe the the- 
oretical framework used to study the interplay of membrane adhesion and do- 
main formation. We first review the behavior of interacting membranes that 
have a homogeneous or uniform composition. We then describe general lattice 
gas models for multicomponent membranes with anchored sticker and repeller 
molecules. Section [3] summarizes the different theoretical methods used to elu- 
cidate the membrane behavior. The remaining sections |4]- [6] represent the main 
part of this review and discuss the interplay of adhesion and domain formation 
for several membrane systems. 
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In section [4j we consider multicomponent membranes with one species of 
sticker molecules and describe several entropic mechanisms that enhance or 
induce the formation of sticker-rich domains within the contact area. These 
mechanisms reflect the large configuration space which the membranes and the 
sticker molecules can explore because of their mobility. In order to elucidate 
these mechanisms, one must distinguish two types of molecular interactions: 
(i) <rans-interactions between the sticker molecules and the second membrane 
or substrate surface; and (ii) cis-interactions between two stickers which are 
anchored to the same membrane. If the cis- interactions between the stickers 
are repulsive and short-ranged, the thermally excited shape fluctuations of the 
membrane induce small clusters of stickers but are not able to initiate the for- 
mation of sticker-rich domains within a sticker-poor membrane matrix. In fact, 
in the latter case, sticker- mediated adhesion occurs only if the sticker concentra- 
tion exceeds a certain threshold value. Sticker-rich and sticker-poor domains axe 
formed if the stickers experience attractive cis-interactions and these attractive 
interactions are effectively enhanced by the thermally excited shape fluctuations 
of the membranes. Furthermore, purely entropic mechanisms for the formation 
of sticker-rich domains arise if the stickers are large compared to the smallest 
wavelength of the bending modes or if the stickers increase the local bending 
rigidity of the membrane. 

Another class of mechanisms for adhesion-induced domain formation is dis- 
cussed in section |5] where we consider the adhesion of multicomponent mem- 
branes which contain both sticker and repeller molecules. If the length of the 
repellers, say Ir, exceeds the length of the stickers, say Is, the size mismatch be- 
tween these two species of molecules favors the formation and growth of sticker- 
rich domains. This can be understood in terms of effective membrane potentials 
with a potential well arising from the stickers and a potential barrier arising from 
the repellers. Phase separation into sticker-rich and sticker-poor (or repeller- 
rich) domains occurs if the potential barrier is sufficiently high. Similar barrier 
mechanisms are also effective if the membranes contains two species of stickers, 
long and short ones. 

Finally, section [6] addresses the time evolution of domains in multicomponent 
membranes which contain (i) stickers and repellers, (ii) short and long stickers, 
or (iii) short stickers, long stickers as well as repellers. In all cases, the effective 
membrane potential exhibits a potential barrier which implies that the initial 
dynamics of the domain formation represents a nucleation process. Thus, in the 
presence of repellers, adhesion is governed by the nucleation of sticker clusters 
or islands. The diffusion and coalescence of these clusters leads to the formation 
of distinct domain patterns at intermediate times. This provides a simple and 
generic mechanism for the observed time evolution within the immunological 
synapse between helper T cells and antigen-presenting cells. 
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2 Modelling of membranes 



In this section, we discuss the theoretical models that will be used in order to 
describe and understand the interplay between membrane adhesion and domain 
formation. This interplay arises from two degrees of freedom: the elastic de- 
formations of the membranes and the spatial patterns of membrane-anchored 
molecules. 

As mentioned in the introduction, we will focus on the contact area of adher- 
ing vesicles or cells and view their unbound segments as reservoirs for the sticker 
and repeller molecules. It is then convenient to use a grand-canonical descrip- 
tion and to describe these reservoirs by chemical potentials for the membrane- 
anchored molecules. In general, the membranes may contain several species of 
such molecules which will be distinguished by the index k with k = 1^ . . . K. 
The membrane concentration of species k is then determined by the chemical 
potential fik- 

Because of the flexibility of the membranes and the lateral mobility of the 
anchored molecules, the membranes can attain many microscopic states within 
the contact area. Each of these states can be characterized by its configurational 
energy or effective Hamiltonian, Tij^ At temperature T, the statistical weight 
of a certain configuration, i.e., the probability to observe this configuration, is 
then proportional to the Boltzmann factor ~ exp(— Ti/fcsT) with Boltzmann 
constant kg. 

The elastic deformations of an adhering membrane will be described by the 
separation field I. The spatial patterns of anchored molecules will be repre- 
sented by the composition variables n. These composition variables are defined 
with respect to an underlying lattice of membrane patches. In this way, the 
well-known theoretical framework of lattice gas models is extended to flexible 
surfaces. [HIKS] 

In general, the membrane-anchored molecules may experience a variety of 
intermolecular forces. We will distinguish between cis- and frans-interactions 
of these molecules. [1] Since two molecules that are anchored to the same 
membrane cannot occupy the same membrane patch, these molecules always 
experience hardcore cis-interactions which are repulsive and short-ranged. In 
addition, these molecules may stick to each another, which corresponds to short- 
ranged attractive cis-interactions, or they may carry electric charges, which 
can lead to long-ranged cis-interactions. In addition, the membrane-anchored 
molecules mediate the trans-interactions between the membrane and the second 
surface. By definition, stickers mediate attractive trans-interactions whereas 
repellers mediate repulsive ones. 

The lattice gas models considered here have two rather useful features: (i) 

■^Here and below, the term 'effective Hamiltonian' is equivalent to the term 'configurational 
energy' which is standard practise in statistical mechanics even though the configurations 
are described in terms of classical, i.e., commuting variables and, thus, do not involve any 
quantum-mechanical degrees of freedom. The Hamiltonians used in this article are 'effective' 
in the sense that they do not describe all molecular details of the biomembranes but focus on 
the relevant degrees of freedom. 
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The hardcore cis-interaction between the anchored molecules, which leads to 
their mutual exclusion within the membrane, is automatically incorporated; 
and (ii) If this hardcore interaction is the dominant cis-interaction, one can 
perform the partial summation over the composition variables in the partition 
function. As a result, one obtains effective membrane models that depend only 
on the separation field 



2.1 Homogeneous or uniform membranes 

2.1.1 Membrane configurations and effective Hamiltonian 

Our theoretical description starts with a homogeneous or uniform membrane in 
contact with another, planar surface with Cartesian coordinates x = {xi,X2)- 
This membrane will be viewed as a thin elastic sheet that exhibits an average 
orientation parallel to this planar surface. If we ignore overhangs, the membrane 
shape can be parametrized by the separation field l{x) which describes the local 
separation of the membrane from the planar reference state with l{x) = 0. 
The elastic deformations of the membrane are governed by two parameters: 

(i) the membrane tension tr, which is conjugate to the total membrane area, and 

(ii) the bending rigidity k that governs the bending energy of the membrane. In 
addition, the membrane experiences direct interactions with the second surface 
that will be described by the interaction energy per unit area or effective mem- 
brane potential KnolO- The membrane configurations are then governed by the 
effective Hamiltonian 

n{i} = n,i{i} + n,n{i} (i) 

which consists of the elastic term 7ioi{^} and the interaction term 7iin{Z}. The 
elastic term has the form 

Wci{^} = j d^x [i a {Vlf + \ K {VHf] (2) 

where the cr-term is proportional to the excess area of the deformed membrane 
and the K-term is proportional to the squared mean curvature of the membrane, 
see Appendix A. The interaction term represents an integral over the effective 
membrane potential as given by 

,{1} = j d?X Knc(0- (3) 

The effective Hamiltonian as given by ([TJ - ([s]) also applies to two inter- 
acting membranes which are characterized by bending rigidities ki and k,2 and 
membrane tensions (Ti and (72, respectively. In this case, the separation field I 
describes the local distance between the two membranes which can both exhibit 
bending deformations and, thus, attain a nonplanar shape, see Appendix A. For 
two interacting membranes, the paramaters a and n in the elastic part ^ now 
represent the effective tension [S] 

a = aii72/{(Ti+ (J2) (4) 



Hi 
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and the effective bending rigidity [IS] 



K = KiK2/(ki + K2). 



(5) 



The interaction potential between the two surfaces always contains a repul- 
sive hardwall potential which ensures that the two surfaces cannot penetrate 
each other and that the separation field I satisfies I > 0. This hardwall poten- 
tial can be implemented in two ways: (i) It may be included into the definition 
of the effective membrane potential via 14nc(0 = °° foi' ^ < (ii) it may be 
embodied by restricting the /-integration in the partition function to positive 
values. In the following, we will use the second implementation and, thus, define 
the partition function via 



2.1.2 Classification of effective membrane potentials 

The effective membrane potential V,no(0 describes the interaction free energy of 
two planar surfaces at uniform separation I = const. In general, the functional 
dependence of V on I will refiect various intermolecular forces such as van der 
Waals, electrostatic, or hydrophobic interactions. [13] In addition, the effective 
membrane potential also depends on external forces or constraints such as an 
applied osmotic pressure or the confinement via another wall that provides an 
upper bound for the separation field I. 

It is useful to distinguish two different classes of membrane potentials cor- 
responding to membrane confinement and membrane adhesion, respectively. 

Membrane confinement is described by potentials 14ne(0 that do not attain 
a finite value in the limit of large I. Simple examples are provided by (i) a 
confining wall at / = L which implies Vmc = 00 for I > L; (ii) a harmonic 
potential as given by ^ G{1 — U)"^; and (iii) an osmotic pressure P which implies 
that V{1) contains the linear term PI for I > 0. 

Membrane adhesion, on the other hand, corresponds to potentials Vme(0 
that (i) have at least one attractive potential well and (ii) attain a finite value 
in the limit of large /. Since we may always replace Vmc by VJnc(0 ~ ^10(00), 
it is sufficient to consider potentials that decay to zero for large I. All effective 
potentials that arise from intermolecular forces have this latter property. 

A confined membrane exhibits critical behavior as one reduces the strength 
of the confining potential. If this potential is symmetric, i.e., if we can define 
a shifted separation field I' = I — lo in such a way that T4io(— = Kno(^')i 
the average value {I') is always zero but the variance {{I' — {l'))^) diverges as 
one reduces the confining potential. This delocalization behavior is obtained 

(i) for a confining wall as one moves this wall to larger values oi I = L and 

(ii) for a harmonic potential as one decreases the strength G of this potential. 
For an asymmetric potential, on the other hand, the confined membrane also 
unbinds from the second surface at Z = as one decreases the strength of the 
confining potential. One example is provided by an external osmotic pressure P: 




(6) 



9 



in the limit of small pressure P, one obtains the complete unbinding behavior 



An adhering membrane, on the other hand, exhibits critical behavior as one 
effectively reduces the attractive part of the effective membrane potential V^e{l). 
Since the membrane configurations are governed by V^e{l)/kBT, this reduction 
can be most easily obtained by raising the temperature T. In the absence of 
any confining potential, the membranes then undergo an unbinding transition 
at a characteristic temperature T = Tu- The nature of this transition depends 
on the functional dependence of the effective membrane potential Vme on the 
separation field The unbinding transition is continuous if Kne(0 has a single 
potential well but no potential barrier. [431 147] Such an effective membrane 
potential arises, e.g., from the interplay of van der Waals and hydration forces. 
If ^mo(0 contains both a potential well and a potential barrier, this transition is 
continuous for sufficiently low potential barriers but discontinous for sufficiently 
high barriers. [Ml US] 

Unbinding transitions have been observed experimentally by Mutz and Hel- 
frich [50 for glycolipid bilayers and by Pozo-Navas et al. |51j for bilayers com- 
posed of two phospholipids. In both cases, the composition of the bilayer mem- 
branes was presumably uniform. 

2.2 Discretized models for uniform membranes 

In order to include the anchored molecules in the theoretical description, we 
will first discretize the uniform membrane. A convenient discretization is pro- 
vided by a square lattice within the planar reference plane, see Fig. [l] The 
corresponding lattice parameter is denoted by a. In this way, the 2-dimensional 
coordinate x is replaced by a discrete set of lattice sites labeled by the index i. 
The membrane configurations are now described in terms of separation fields U 
associated with the lattice sites i, and the membrane is divided up into discrete 
membrane patches, each of which has projected area a?, see Fig.[l] 

Since the elastic part of the effective Hamiltonian depends on the derivatives 
of I with respect to the coordinates Xi and X2, we have to discretize these deriva- 
tives as well. For the excess area term ^ (VZ)^, we will use the discretization 

{VdkT = [K^i + a,X2) - l{xi,X2)f + [l{xi,X2 + a) - l{xi,X2)f (7) 

where xi and X2 denote the Cartesian coordinates of the lattice site i. Likewise, 
the discrete Laplacian is taken to be [47] 

V^k = l{xi +a, X2) + l{xi — a, X2) + l{xi,X2 + a) + l{xi,X2 — a) —Al{xi,X2) (8) 
The elastic Hamiltonian now has the form 



(9) 



and the interaction Hamiltonian becomes 




(10) 
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which represents a summation over all membrane patches i with potential en- 
ergies 

Vik) = a^V^,{l{xi,X2)) (11) 

where xi and X2 denote again the Cartesian coordinates of the lattice site i. 
Note that V{li) and Vnio(0 have the dimensions of energy and energy per area, 
respectively. 

Realistic estimates of the entropy and free energy of the membranes require 
that the lattice constant a is equivalent to the smallest possible wavelength for 
bending fluctuations of the membranes. Computer simulations with molecular 
membrane models indicate that this size is somewhat larger than the thickness 
of the lipid bilayer and of the order of 6 nm 113] ■ 

2.3 Lattice gas models: General form 

Next, we include the membrane- anchored molecules, that may act as stickers 
or repellers, into the theoretical modelling. In general, the membrane may 
contain K different types of such molecules which will be distinguished by the 
index k = 1,...^K. Since all membrane-anchored molecules undergo lateral 
diffusion along the membrane, these molecules can form many different spatial 
patterns. In order to describe these patterns, we now introduce composition 
variables rii for all lattice sites i. Each composition variable can attain the values 
rii — 0,1, ...,K. If the membrane patch i contains the membrane-anchored 
molecule of type k, this patch is characterized by = fc. If the patch does not 
contain any of the K membrane-anchored components, the composition variable 
rii has the special value = 0. 

In the absence of the second membrane or surface, the concentrations of 
the K species of membrane-anchored molecules are governed by K chemical 
potentials Hk with k = 1,2, . . . , K. In addition, the cis-interaction between one 
molecule of species k located at lattice site i and another molecule of species k' 
located at site j is described by the pair-potential W^j'' which is negative for 
attractive cis-interactions. Thus, the configurations of the composition variables 
rii are governed by the cis-interaction part of the effective Hamiltonian as given 

by 

K K K 

^ciM = - E E ^fc'^fc^". + E E E 4,„,4',„,M/^''' (12) 

i k=l (ij) k=l k' = l 

with the Kronecker symbol 6k^n which is defined by Sk^n = 1 for k ~ n and 5k,n = 
otherwise. The symbol {ij) indicates a summation over all pairs of lattice 
sites i and j. Note that the chemical potential term alone already embodies the 
hardcore interactions between two neighboring membrane-anchored molecules 
because of the underlying lattice. 

The configuration of the adhering membrane is now described both by its 
separation field U and by its composition variables rii. These two degrees of 
freedom are governed by the effective Hamiltonian 

n{l,n} = H^x{l,n} + nin{l,n} (13) 
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where the elastic part is now given by 

7iei{/, - a {^dkf + 2^ {^dkf] (14) 

i 

where the bending rigidities Ki depends, in general, on the composition variable 
rii. In the following sections, we will mostly consider the simplified situation 
characterized by an ^-independent bending rigidity Ki — k; an exception is 



section 4.5 where the presence of sticker molecules leads to more rigid membrane 
patches. 



The interaction part of the effective Hamiltonian (13) has the more general 
form 

nin{l,n} = nci{n} + HtAlri} (15) 



which consists of the cis-interaction part Tid as given by ( 12 ) and the trans- 
interaction part Titr which describes the various trans-interactions between the 
two surfaces as mediated by the different membrane-anchored molecules. 

The membrane patch i experiences the trans-interaction Vkili) if the patch 
contains a membrane-anchored molecule of type k, and the trans-interaction 
Vo(^i) if this patch does not contain any membrane-anchored molecule. The 
total trans-interaction part of the effective Hamiltonian is then given by 

K 

Wtr{Z,ri} = ^^4,„.^fe(Z,) (16) 

i fe=0 

with the Kronecker symbol (5fc.„ as before. 

So far, we have not specified the projected area of the membrane-anchored 
molecules. If this projected area does not exceed the area a? of the membrane 
patches, the short-ranged repulsive cis-interactions between the membrane-anchored 
molecules are incorporated by the underlying lattice and the composition vari- 
ables rii. Thus, if the membrane- anchored molecules do not experience at- 
tractive or long-ranged repulsive cis-interactions, the interaction Hamiltonian 
attains the simple form 

K 

7ii„{?,n}-^^4,„jyfc(/,)-^fe] with ^0 = 0. (17) 

i k=n 

This form is particularly useful since one may now perform the partial sum- 
mation over the composition variables in the partition function as will be 
explained in more detail below. Furthermore, we will typically assume that 
the trans-interaction Vo{li) for membrane patches without any stickers or re- 
pellers corresponds to a short-ranged repulsion that can be incorporated into 
the hard wall repulsion that ensures > for all i. 

2.4 Membranes with sticker molecules 

The simplest example of adhesion mediated by membrane-anchored molecules 
is provided by membranes that contain only one species of sticker molecules. 
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Figure 1: Segment of deformable membrane in contact with a homogeneous and 
planar surface (such as another membrane supported on a planar substrate). 
The planar surface provides a 'reference' plane which is discretized into a square 
lattice. The deformable membrane is composed of two types of patches: Mem- 
brane patches without a sticker molecule are grey wheras membrane patches 
with such a molecule are black. The configurations of the deformable membrane 
are described (i) by the separation field k, which measures the local separation 
of the two surfaces, and (ii) by the composition variable n,j with n,j = for 
sticker-free patches and = 1 for patches with a sticker molecule. 



[TJ[51[3] In this case, the composition variables attain only two values: — 1 
corresponding to a sticker in membrane patch i, and ti^ = corresonding to no 
sticker in this patch, see Fig.[l] Therefore, the composition variables rii are now 
equivalent to occupation numbers for the sticker molecules. 

Inspection of Fig.[l]shows that the membrane models considered here repre- 
sent effectively two-component systems. It is important to note, however, that 
the two components do not correspond to individual lipid or sticker molecules, 
but to lipid bilayer patches of size with or without a sticker molecule. In the 
absence of the stickers, the lipid bilayer patches need not consist of a single lipid 
species but are assumed to have a homogeneous or uniform composition. In this 
sense, the models described here correspond to multi-component membranes 
that can have many nonadhesive components but only one adhesive component. 

In section [4] we will consider several models of multi-component membranes 
with one species of membrane-anchored stickers that interact with an effectively 
homogeneous or uniform membrane. Such a situation arises if the stickers of 
the multi-component membrane bind directly to a lipid bilayer with uniform 
composition, or if the stickers bind to ligands in the second membrane that are 
present at high concentrations. These models can be extended to more complex 
situations in which both membranes contain sticker molecules ^ [3] . 

The effective Hamiltonian for a multi-component membrane that interacts 
with a second, homogeneous membrane via one species of stickers can be ob- 
tained from the general expression ( |13[ ) with K ~ 1. One then obtains the 
effective Hamiltonian 7i{Z, n} = Hc\{l, n} + Ti.ci{n} + 7itr{^, n} with the elastic 
part Hci{l,n} as given by p4| . 

The cis-interaction part of the effective Hamiltonian now has the relatively 
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simple form [H [21 E] 



Hci{n} ^ - ^^J■ni + ^ W^jTiin^ (18) 

where {ij) indicates a summation over all pairs of lattice sites i and j, and the 
parameters /i = /ii and Wij = W^j^ represent the sticker chemical potential 
and the sticker-sticker pair potential, respectively. Attractive and repulsive 
pair-potentials are described by Wij < and Wij > 0, respectively. 

The trans-interaction part of the effective Hamiltonian is now given by [1] 

HtAl, nj^Yl - ^^)Vo{k) + n,V,{k)] , (19) 

i 

where Vq is the trans-interaction of a membrane patch without sticker whereas 
Vi is the trans-interaction of a patch with sticker. 

The interaction part of the effective Hamiltonian again simplifies if the pro- 
jected area of the sticker molecules does not exceed the patch area and if two 
sticker molecules do not experience attractive cis-interactions or long-ranged 
repulsive ones. If one can also neglect the trans-interaction Vq for a sticker-free 
membrane patch, the interaction Hamiltonian is simply given by |i2j [3] 

Hi„{;,n} = ^n,[yi(z,)-/i]- (20) 



In section 4.4 we will also consider sticker molecules which cover several 
membrane patches. This implies that the projected area of the sticker molecule 
exceeds the patch area and that the smallest separation of two sticker molecules 
is larger than the lattice parameter a. This situation will be described by repul- 
sive pair-potentials Wij > which extend to next-nearest neighbor lattice sites 
in order to prevent the overlap of adjacent sticker molecules. 



2.5 Two types of membrane-anchored molecules 

Another relatively simple case is provided by multi-component membranes with 
two species of membrane-anchored molecules ^31 13 E] • In section [s] we will 
discuss membranes with one type of sticker and one type of repeller molecule. 
Likewise, one might consider membranes with two species of sticker molecules, 
say short and long ones. In these cases, the composition variables rii can attain 
three different values = 0, 1, and 2 where the value = corresponds to 
a membrane patch without any membrane-anchored molecule as before. For 
stickers and repellers, the values Ui = 1 and rii = 2 are taken to indicate a 
sticker and a repeller molecule, respectively. 

It is again useful to consider the simplified situation in which (i) the hardcore 
interaction between the stickers and repellers represents their dominant cis- 
interaction, and (ii) the trans-interaction Vo(/) for membrane patches without 
a sticker or repeller molecule is well described by the hardwall potential which 



14 



ensures I > 0. The interaction part of the effective Hamiltonian then has the 
simple form O [S] 



Hi^il, n} = Yl (^-(^») - + ^2.«. (Vrih) - Mr)] (21) 



as follows from (17 1 with the sticker chemical potential fig = /ii, the repeller 
chemical potential fir = M2, the trans-interaction Vs = Vi mediated by the 
stickers, and the trans-interaction Vr = V2 mediated by the repellers. 

The equilibrium phase behavior for a membrane with anchored stickers and 
repellers will be reviewed in section |5.1[ The adhesion dynamics of a multi- 
component vesicle with stickers and repellers is the topic of section |6.1[ In 
the latter case, the unbound membrane segment which acts as a reservoir for 
the sticker and repeller molecules will be taken into account explicitly, and the 
total number of membrane-anchored molecules within the bound and unbound 
membrane segment will be kept constant. 

As mentioned, the linear size a of the membrane patches corresponds to the 
smallest possible wavelength, or 'cut-off' length, for bending fluctuations of the 
membrane. This length is somewhat larger than the bilayer thickness [55J [53] 
and affects the entropic, fluctuation-induced mechanisms for phase separation 
discussed in section [4] However, the barrier mechanisms for the lateral phase 
separation of membranes with stickers and repellers as discussed in section |5] 
are rather insensitive to the precise choice of the cut-off length a, as follows 
from scaling estimates [5] . We will use this property in order to choose a larger 



patch size for the T cell adhesion model [5] as reviewed in subsection 6.2 This 
choice is convenient since one has to consider relatively large cell-cell contact 
areas, and two different species of stickers which have a rather different linear 
size. This larger choice for the size of the membrane patches also implies that 
these patches can contain more than one sticker or repeller. 



3 Theoretical methods 

3.1 Monte Carlo simulations 

3.1.1 Simple sampling and importance sampling 

In classical statistical mechanics, systems are characterized by their configura- 
tions {y} and their configurational energy, or Hamiltonian, Ti.{y}. The thermo- 
dynamic properties of the system can be expressed as averages 

(A(y)) = l^A(y)e-«{^>/'=-^ (22) 
y 

where Z is the partition function 

Z = ^e-«{^>/'=«^ (23) 
y 
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For many systems of interest, partition functions and thermodynamic averages 
cannot be calculated exactly. These systems have to be studied with approxi- 
mate analytical methods or with numerical algorithms. Numerical methods face 
the problem that the systems of interest have usually many degrees of freedom 
and, hence, a large configurational space. These methods therefore can only 
probe a subset or 'sample' of the full configurational space. 

In the simple or random sampling Monte Carlo method, a subset {yi, . . . , um} 
of the configurations space is randomly selected. The average of a certain phys- 
ical quantity, say A, is then estimated by 



m = ^':^r"L^..,. ^ (24) 



In practice, the random sampling method works only for systems with a flat 
energy landscape for which all configurations have the same energy; one example 
is provided by random walks on a lattice |54j . Simple sampling is ineffective 
for other systems where only a relatively small fraction of configurations yi has 
a low energy TL{yi} and, hence, a large Boltzmann weight exp[— Tilj/ij/ZcsT]. 
These configurations then dominate averaged quantities {A{y)) as in ( [22| , but 
are hard to find from random sampling. 

Another sampling method, which is more efficient than random sampling, 
is the importance sampling Monte Carlo method. With this method, config- 
urations are not randomly selected but 'generated' successively via a Markov 
process. One important constraint on this Markov process is that the system 
relaxes, for long times, towards the equilibrium distribution 

= ^e-^^y^y^-^ (25) 

In this Markov process, a configuration is generated from the preceding config- 
uration yi-i with a certain transition probability W{yi-i yi). The transition 
probabilities have to fulfill the detailed balance conditions 

Peq{yi) W{y, ^ = W{y, yi) (26) 

which ensure that the distribution of generated states relaxes towards the equi- 
hbrium distribution [511 ES] . 

Using the expression ( |25| for the equilibrium distribution, the detailed bal- 
ance condition can be rewritten in the form W{yi yj)/W{yj yi) = 
exp[— (7i{yj} — Hlyi}) /kBT] . This means that the ratio of the forward and 
backward transition rates between the states i and j only depends on the energy 
difference Hj — Hi- One choice for the transition rates that obeys detailed bal- 
ance is provided by the Metropolis algorithm which is defined by the transition 
probabilities |54| 

Wiv^^y,) = e-(«te>-«{^'})/'=«^ forH{y,}-7^{2/.}>0 
= 1 otherwise . 
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3.1.2 Membrane simulations 



In order to illustrate how Monte Carlo methods can be used to study the discrete 
membrane models introduced in subsection |2.2| above, let us first consider two 
homogeneous and tensionless membranes with interaction potential V{1). It is 
convenient to introduce the rescaled separation field 



z = {lla)^K/kBT (28) 

where a is the lattice parameter (or linear patch size) of the discretized mem- 
brane and K is the effective bending rigidity. Introducing the rescaled separation 
field z simplifies the notation and reduces the number of independent param- 
eters. In terms of the rescaled separation field z, the Hamiltonian for the two 
membranes has the simple forn:]^ 

^{^} = E^{(^'^^»)' + ^(^»)} (29) 

■i 

where is the discretized Laplacian. If «1, i2, i3 and i4 denote the four nearest 
neighbors of site ? on a square lattice, see Fig. [2] the discretized Laplacian can 
be written as A^z^ = zn + zii -\- za -\- zn — Azi. Usually, a 'new' membrane 
configurations in the Markov process is generated from an 'old' configuration by 
attempting a local move Zi — Zoid — * -^i = -Zncw for a randomly selected lattice 
site i. Choosing a new value for the rescaled membrane separation z^ at site i 
affects not only the discretized Laplacian at this site, but also at the four nearest 
neighbor sites il, i2, z3 and lA. Therefore, the energy difference between the 
'old' and 'new' membrane configuration is given by 

AH = Hncvj — Hold = '^|zi=2„o„ — '^|zi=2old 

= i [{Adz,f + {Adzaf + (Arfz,2)2 + {Adz.^f + (Adz.i)^] - 
\ [(Arfz,)2 + (Adza)^ + (Arfz,2)2 + {Aaz^f + {Adz,^f\ ^ 

^^(^ncw) - V^(^old) 

(30) 

If i5 to i8 denotes the next-nearest neighbors of site i, and i9 to il2 the 
next-next-nearest neighbors (see Fig. [2]), the energy difference can be written as 

'Wnow — Wold = 10 (Zncw ^ ^old) + [^now " 2old][— 8(zil + Zi2 + ^iS + ^ia) 

+ 2{zi5 + Zid + za + Zis) + Zi-s, + z^o + ^^i + zi-y-^ -I- y(z„ow) - ^^(^^oid) (31) 



In the Metropolis dynamics ( 27 1 , the move is always accepted if Tinow — "^oid is 
negative, but is only accepted with probability exp[— (Tincw — "^oid)] if "^now — 
Tioid is positive. Of course, these local moves have to be attempted at all sites or 
membrane patches, and typically a large number of Monte Carlo steps per site is 
required to reach the equilibrium distribution from a given initial configuration. 



*Here and below, H{2:} and V{zi) are dimensionless quantities given in units of the thermal 
energy ksT. 
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Figure 2: Neighbors of a site, or membrane patch, located in the center of the 
figure. Here, 1 to 4 are the nearest neighbors, 5 to 8 the next-nearest neighbors, 
and 9 to 12 the next-next-nearest neighbors. 



It is important to choose new values ^;new for the rescaled separation in an 
unbiased way, e.g., using the rule 

Znew = Zold + SzC (32) 

where Sz is the step size, and C is a random number between -1 and 1. Steps with 
negative values for z,icw art; rejected since the two membranes cannot penetrate 
each other. 

The relaxation time in units of Monte Carlo steps depends on the step size 
6z. If 5z is large, then only a small fraction of Monte Carlo steps will be 
accepted and the relaxation is slow. If 6z is small, most Monte Carlo steps will 
be accepted, but the relaxation is also slow because it takes a large number of 
steps per site to obtain a significantly different configuration. Usually, a suitable 
value for 6z can be found by trying to obtain an average acceptance rate between 
0.4 and 0.5 for the Monte Carlo moves. 

The simulation of multicomponent membranes requires Monte Carlo moves 
of the concentration field n in addition to the moves of the separation field 
z. In the grandcanonical ensemble, for example, a simple Hamiltonian for a 
membrane with stickers is given by 

n{z, n} = Heiiz} + [V{zi) - m] (33) 

i 

where /x is the chemical potential of the stickers. The local concentration rii 
adopts the values or 1 indicating the absence or presence of a sticker at site i. 
The Monte Carlo move rij = rioia — > ?^new = l — rii attempts to remove stickers 
from sites with rij = 1, and to add stickers at 'empty' sites with rij = 0. The 
energy difference for these moves is 

Wnew - Wold = - = (1 - 2n,)[V{z,) - fi] (34) 
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3.2 Free energies of adhesion for homogeneous membranes 

Thermally excited shape fluctuations lead to an entropic repulsion between 
membranes |56j . If the membranes are bound together by attractive trans- 
interactions, thermal fluctuations can cause an unbinding transition at a criti- 
cal temperature At temperatures T below T^, the attractive interaction 
between the membranes dominates, and the membranes are bound. At tem- 
peratures above Tc, the entropic repulsion dominates, and the membranes are 
unbound. Renormalization group calculations show that the unbinding tran- 
sition is continuous, i.e., the average membrane separation (l) \T — T^j"^ 
diverges continuously when the critical temperature is approached from below . 
This continuous behavior has been confirmed by Monte Carlo simulations [17] . 

Monte Carlo simulations can also be used to determine the free energy of 
adhesion [21 [3] . The discretized Hamiltonian for two homogeneous membranes 



has the form as given by (|9| and (10 1. In the following, the lateral tension a of 
the membranes is assumed to be zero and the attractive interaction potential of 
the two membranes is taken to be a square well potential with potential depth 
ksTu < and potential range ly as given by 

V{k) = kBTue{U-k) =kBTu for < /i < 

= for k > U. (35) 



After introducing the rescaled separation field z = [1/ a)yj n/kBT , the system 
is described by two dimensionless parameters, the potential depth u and the 
rescaled potential range = (l^j a)^ nlksT ■ The hard wall potential which 
restricts the separation field li to positive values is incorporated via a lower 
bound for the /^-integration. 

In general, it is difficult to determine partition functions and, thus, free 
energies via Monte Carlo simulations. However, the free energy of two adher- 
ing membranes can be determined via the standard method of thermodynamic 
integration. The quantity of interest here is the contact probability 

Pfc = {B{U ~ k)) (36) 

which represents the average value for the fraction of bound membrane seg- 
ments, i.e., membrane segments with local separations < 1^ < ly. The contact 
probability is related to the free energy density 



-ksTln 



/>oo 

n/ di 



^H{l}/kBT^^^ (37) 



which represents the free energy per membrane area A, via 

kB-i ou 

The contact probability can be determined by Monte Carlo simulations, see 
Fig. [3j For a given rescaled potential range z„, the free energy of adhesion is 
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Figure 3: Contact probability Ph of a homogeneous, tensionless membrane as a 
function of the depth |m| of the square-well potential for the rescaled potential 
ranges z.^ — 1, 0.5, 0.3, and 0.1 (from left to right). 



then obtained from the Monte Carlo data for the contact probability by the 
integration 

F - i^„6 - ^ du'Pfc(u') . (39) 

Here, F^h denotes the free energy density of the unbound membrane, and u„ 
is the critical potential depth. Close to the critical point with P5 = 0, the 
correlation length and relaxation time of the membranes diverge. Therefore, 
small values of Pf, cannot be determined reliable in the simulations. But since 
P\y is proportional to (u — m.^ ) close to the critical potential depth , the critical 
points and the full functions Pf,(u) can be obtained using a linear extrapolation 
of the data shown in Fig. |3] 
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3.3 Effective potential in the absence of cis-interactions 



As explained in section |2.4[ the simplest Hamiltonian for a discrete membrane 
with stickers has the form 



H{Z,n}=?^,i{?} + ^n,;[F(Z,)-/x] 



(40) 



where the elastic part Tioi is given by ( 14 1 with uniform bending rigidity k. 



The stickers have the size of the lattice parameter a and cis-interactions between 
the stickers are absent. Since the effective Hamiltonian in (40 1 is linear in n. 



the partial summation over the sticker degrees of freedom can be performed in 
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Figure 4: Unbinding lines for the rescaled potential ranges = 1, 0.5, 0.3, and 
0.1 (bottom to top). The membranes are bound at large sticker concentrations 
X and binding energies \U\, and unbound at small sticker concentrations and 
binding energies. 



the partition function 



n 



dk 



n E 



ni=0,l. 



-H{l,n}/kBT 



(41) 



which leads to 
Z = 



dh 



[n^i{i}+T,iV.i(h))/kBT 



(42) 



Here, is the number of lattice sites. For a square-well sticker potential V{li) = 
U6{lv — h), the effective potential has the form 

1 + Jlx-U]/kBT 

Vefik) = -ksTln -Z___0Q^ _ = u,f9{l, - k) (43) 

This effective potential is again a square-well potential. It has the same poten- 
tial range as the sticker potential, and an effective potential depth f/g/ which 
depends on the chemical potential n and the binding energy U of the stickers. 

Summing out the stickers degrees of freedom in the partition function thus 
leads to the problem of homogeneous membranes with square-well interaction. 
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The free energy per unit area T = —{ksT/A) InZ therefore can be obtained as 



in the previous section 3.2 
keT 



T - 



In 



1 



du'Pbiu') 



(44) 



Here, u — Uef/ksT is the dimensionless potential depth u, and Pb{u) is the 
contact probabiHty for homogeneous membranes. The effective potential depth 
Uef is defined in eq. (43 1. The sticker concentration X = —{dF/djj) = {ni)/a^ 
then follows as 



X = 



(1-n) 



1 -|- QtJ./kB' 



^[tJ.-U]/kBT 



(45) 



To obtain the critical sticker concentration X^, at which the membranes unbind, 
we have to determine the critical chemical potential /i* of the stickers. This 
critical chemical potential can be derived from the equation m» = Uef{lJ^*) /ksT: 



knT 



= ln 



,\u,\ 



p\U\/kBT _ g|«.| 



(46) 



At the continuous unbinding point, the contact probability Pb is zero. The 
critical sticker concentration therefore is 



X:^ 



^t^./kBT 



1 



a2(i + eA'-/feB'J^) a^{e\U\/kBT 



(47) 



The resulting unbinding lines for various rescaled range potential Zy are shown 
in Fig.|4] The lines follow from eq. (47l and the critical values m* = — 1.25±0.05 
for z„ = 0.1, = -0.47 ± 0.03 for z„ = 0.3, = -0.24 ± 0.02 for z„ = 0.5, 
and = —0.075 ± 0.010 for — 1. These values are obtained from linear 
extrapolation of the MC data for the contact probablility Pb shown in Fig. [3] 
The two-component membrane considered in this section does not exhibit 



lateral phase separation. This is a direct consequence of eq. (45 1 and the conti- 
nuity of the contact probability Pb- 

3.4 Variational (mean-field) theory for cis-interactions 

As mentioned before, a sticker molecule might also experience cis-interactions 
with other stickers in the same membrane. In general, these cis-interaction 
can be attractive or repulsive and may be short- or long-ranged. If the trans- 
interaction Vq for membrane patches without stickers are again purely repulsive 
and short-ranged, the effective Hamiltonian has the form 



H{1, n} = -Heiil} + ^ [V{k) - A*] 



(48) 



as follows from (^18| and (19). Here, Wij is the cis-interaction energy between 



two stickers in the membrane patches i and j. The summation index (ij) in- 
dicates a summation over all pairs of membrane patches. In the following, we 
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consider the short-range interaction energy 



Wij = W for nearest neighbors i, j 

= otherwise (49) 

which only has a nonzero value W for two stickers in adjacent membrane patches. 
The cis-interaction term in the Hamiltonian (58 1 contains an expression Uirij 



which is quadratic in the concentration field. Therefore, a direct summation of 
the sticker degrees of freedom as in the previous section is not possible. Ap- 
proximate analytical methods addressing such situations are mean-field theories, 
which are sometimes also called self-consistent field theories. Mean-field theo- 
ries can be derived in a systematic way from the variational principle |57j 

^ <^o + j{n-no)o (50) 

Here, Tio typically is a Hamiltonian which is linear in the considered field. For 
such a Hamiltonian, the degrees of freedom can be summed out exactly in the 
partition function. The variational principle states that the free energy J- of 
a system with Hamiltonian Ti. is smaller than or equal to the free energy of a 
system with Hamiltonian Tio plus the average value of the 'perturbation' TL—TLo- 
This average value is calculated in the system with Hamiltonian TLo, which is 
indicated by the subscript of the brackets (• • • )o- In the present situation, the 
Hamiltonian Ho is 

Uoil, n} = Hail} + n^{V{h) - fi + B) (51) 

i 

where i? is a variational parameter and A the membrane area. The effective 
Hamiltonian Tio is linear in the sticker concentration field n and corresponds 
to stickers with W — and the shifted chemical potential fi — B. Therefore, to 



obtain the free energy Fo we simply have to replace by /i — i? in eq. (44 1. 
Evaluation of {TC — Ti.o)o now leads to 

^ < •^o-^^(n.)o-h^^l^„(n,n,)„ (52) 

< J^o- ^{n^)o+'^{n^)l = J'> (53) 

In deriving the second inequality, we made use of the relation ( 
The relation results from the fact the correlation function {ninj)o— {ni)^ is non- 
negative because of fluctuation-induced attractive interactions between bound 
stickers. These fluctuation-induced interactions will be discussed in section 33] 
The factor 2 in the term 2W/a'^ is the number of nearest-neighbor patches mul- 
tiplied with 1/2 to avoid a double-count of neighbor pairs. 

The variational principle ( |50[ ) tells us that the best approximation for the 
free energy J- is the minimum of J-^ with respect to the variational parameter 
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B. Thus, we are looking for a value of B satisfying 



^ = ^ - in.)o + {mn.)o -Bf-^ = (54) 

Since dTo/dB = —dJ-o/dji = {ni)o, the minimization leads to the self-consistency 
equation 

B = AW{n,)o (55) 

Another name for this equation is 'mean-field equation' since it is also obtained 
in the simple approximation jWijUiUj « jWijni{nj) — 4W{nj)J2i''T-i 
in which the effect of the four neighbors rij on rii is taken into account as the 
'mean field' 4:W{nj). However, to derive the phase behavior, we also need the 



free energy ( 53 1 . A first-order phase transition corresponds to two solutions 



i?i(/i) and -02 (m) of the self-consistency equation (551. At the transition point, 
the free energies ^>(i3i(/i)) and J^>(_B2(/i)) of the two coexisting phases have 
to be equal. From the latter condition, one can derive the transition value fit 
of the chemical potential. 



4 Entropic mechanisms for domain formation 

4.1 Entropic interactions between bound stickers 

Membranes that are confined by external forces or constraints loose configura- 
tional entropy. This entropy loss is proportional to the projected area of the 
membrane and may be viewed as as effective repulsive force as first proposed by 
Helfrich [56l|58]. For a fiuctuating membrane in contact with another surface, 
this fluctuation-induced repulsion behaves as 

Vfi^CfiiksTf/Kp (56) 

for large average separations I between the membrane and the second surface 
where k is, in general, the effective bending rigidity for the separation field I 
of the two surfaces. The dimensionless coefficient c// has been determined by 
Monte Carlo simulations and is found to be of the order of 0.1 j47l \59\ |60l l6T| . 
The precise value of Cfi depends on the confinement of the membranes and 
varies from c/; ~ 0.08 for a single membrane between parallel walls |59j to 
Cfi = 0.115 ± 0.005 for a membrane that is pushed by an external pressure 
against a wall [47l IBTl [1]. In a more systematic treatment, this entropically 
induced interaction arises from the renormalization of the hardwall potential 
Vhw(0 with yhw(0 = for / > and Fhw(0 = oo for / < [43,. 



The fiuctuation-induced repulsion as given by ( 56 ) has also been used to esti- 
mate the shape of a fiuctuating membrane that is pinned to another membrane 
or surface by some external constraint ESI EH Ei] • A detailed comparison 
between self-consistent calculations based on an effective Hamiltonian Ti.{l} with 
the membrane potential Vmc(0 = Cfi{kBT)'^ /kP and Monte-Carlo simulations 
of the ffuctuating membrane subject to local pinning forces has shown that the 
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Figure 5: (Left conformation) A fluctuating membrane is pinned to a surface 
at two different spots by stickers that are far apart from each other. (Right 
conformation) If the two stickers are close together, the membrane is effectively 
pinned only at one spot, which costs less conformational entropy of the fluctuat- 
ing membrane. The membrane fluctuations thus induce an attractive interaction 
between the stickers. 



dimensionless coeflicient Cfi depends both on the boundary conditions imposed 
by the pinning forces and on the physical quantity under consideration |65j . 

The membrane fluctuations also induce interactions between bound stick- 
ers. These interactions are attractive, as one can intuitively see by inspection 
of Fig. |5] Suppose we have two fluctuating membranes with an average sepa- 
ration much larger than the sticker binding range. Locally clamping the two 
membranes together by a single bound sticker then costs a certain amount of 
work w against the entropic repulsion of membranes. If we clamp the mem- 
branes together by two stickers which are far away from each other, we have to 
compress the membranes at two spots, which costs 2w. But if the two stickers 
are very close together, we only have to compress the membranes at a single 
spot, with an entropic cost around w. Therefore, the right configuration in 
Fig. [5] is entropically preferred, since the left configuration constrains, or sup- 
presses, the membrane fluctuations more strongly. The fluctuations thus induce 
an attractive interaction between the bound stickers. 

In principle, the strength of these attractive interactions between bound 
stickers can be determined by integrating out the membrane fluctuations. More 
precisely, one has to integrate out the degrees of freedom for the separation fleld 
/ in the partition function, for fixed locations of the bound stickers. In practice, 
integrating out the membrane fluctuations is not possible in a rigorous way, 
at least for an arbitrary sticker configuration. Approximate scaling arguments 
have been used to estimate the entropic interactions for a regular array [16] and 
an isolated pair of bound stickers |66j . Since many-sticker interactions such as 
screening effects clearly are important, the fluctuation-induced interactions can 
again not be obtained as a simple sum of pairwise interactions. 

In the following, we will focus on approaches which directly address the 
phase behavior of membranes with stickers. The central question in this section 
will be if the fluctuation-induced interactions between the stickers can be strong 
enough to induce lateral phase separation. We will see that the answer to this 
question depends on the type of stickers. 
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Figure 6: Typical Monte Carlo configurations of a two-component membrane 
with small, flexible stickers in the absence of cis-interactions. Bound stickers 
are black, unbound stickers gray, and membrane patches without stickers are 
white. The average separation of the membranes increases continuously as the 
sticker concentration is reduced. The bound stickers tend to form small clusters 
because of fluctuation-induced interactions. However, the membrane does not 
phase separate. 
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4.2 Small flexible stickers without cis-interactions 



In the simplest case, the stickers are small and flexible and have no cis-interactions. 
In the grand-canonical ensemble, the Hamiltonian then has the form |TJ |3] 

n{l,n} = nei{l} + J2n,[Vik)-n] (57) 

i 

where Ti.ei{l} is the elastic energy of the membranes, and V(li) and fj, are the 
trans-interaction potential and the chemical potential of the stickers. In section 
|3.3[ we saw that the sticker degrees of freedom in the partition function can be 
summed out exactly, because the Hamiltonian is linear in the concentration field 
n. Summing out the sticker degrees of freedom leads to a partition function for 
homogeneous membranes interacting with an effective potential Vef{li). If the 
sticker potential V{li) is a square-well potential, then Vefih) is a square- well 
potential, too. It has the same potential range and an effective potential depth 
Uef which depends on the chemical potential fi and the binding energy U of the 



stickers, see eq. ( 75 1 



The unbinding transition of homogeneous, tensionless membranes with at- 



tractive and short-ranged interactions is continuous |43], see also section 3.2 
This implies that the contact probability Pf, , the fraction of bound membrane 
segments, continuously goes to the zero as the unbinding point is approached. 



The sticker concentration X given in eq. (45 1 then is a continuous function of 
the chemical potential /i of the stickers. Therefore, lateral phase separation 
in sticker-rich and sticker-poor domains does not occur |5] [3]. In the grand- 
canonical ensemble, lateral phase separation corresponds to a discontinuity or 
jump in X(/i) at a certain value of the chemical potential. 

The Monte Carlo snapshots shown in Fig. [6] illustrate the continuous un- 
binding of the membranes. The average separation of the membranes increases 
continuously with decreasing sticker concentration, since membrane fluctuations 
become more and more pronounced. The membrane fluctuations also lead to 
small clusters of bound stickers. However, the fluctuation-induced interactions 
between the stickers do not cause lateral phase separation. 

4.3 Stickers with cis-interactions 

If the stickers now interact via cis-interactions, the Hamiltonian has the form 

m 

H{1, n} - Heiil} + 51 i^ih) -tA + Y. "^"J (58) 

In the following, we consider the attractive and short-ranged cis-interactions 

Wij — W for nearest neighbors i, j 

— otherwise (59) 

with characteristic interaction strength W < 0. 
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Figure 7: Mean-field phase diagrams for small, flexible stickers with attractive 
nearest-neighbor cis-interactions. Here, W is the strength of the cis-interactions, 
U is the binding energy, and X is the concentration of the stickers. The phase 
diagrams at W/ksT — —0.33, —0.4, and —0.6 contain shaded coexistence re- 
gions in which the membrane phase separates into a bound sticker-rich and an 
unbound sticker-poor phase. 



Clearly, large absolute values \W\ of the interaction strength will lead to 
lateral phase separation into sticker-rich and sticker-poor domains. The cen- 
tral question here is how membrane fluctuations affect the critical interac- 
tion strength for lateral phase separation. An important reference value is 
the critical interaction strength Wc for 'non-fluctuating' membranes with con- 
stant separation li. In the absence of membrane fluctuations, phase separa- 
tion can only occur for interaction strengths with \W\ > \Wc\. For constant 
membrane separation li, the Hamiltonian (58 1 reduces to the Hamiltonian of 
a 2-dimensional (2d) lattice gas. This lattice gas has the critical interaction 
strength Wc/hgT = — 21n(l -I- ^/2). Mean-field theories, however, systemati- 
cally underestimate the critical interaction strength and lead to Wc/ksT — —1. 

Mean-field phase diagrams for fiuctuating membranes are shown in Fig. [7] 
The rescaled potential range of the stickers here is = 0.1. The mean-field 
theory for the concentration field n is described in detail in section |3.4[ The 
phase diagram at the interaction strength W/kgT = —0.3 contains a bound and 
an unbound phase which are separated by a single line of continuous unbinding 
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Figure 8: Comparison of Monte Carlo (data points) and mean-field phase di- 
agrams (lines and shaded coexistence region). The cis-interaction strength is 
W = 0.5Wc, and the rescaled potential range is = 0.1. Here, Wc is the criti- 
cal interaction strength of the 2d lattice gas, which has the value Wc/ksT = — 1 
in mean-field theory and Wc/ksT = — 21n(l -I- \/2) for the Monte Carlo simu- 
lations. 



transitions. The membrane is bound for high concentrations or high binding 
energies \U\ of the stickers, and unbound for low concentrations X or binding 
energies. For larger absolute values \W\ of the cis-interaction strength, the dia- 
grams contain two-phase regions, shaded in gray. In the two-phase regions, an 
unbound phase with a low concentration of stickers and a bound phase with a 
higher sticker concentration coexist. The coexistence regions end in tricritical 
points. For binding energies \U\ below the tricritical value \Utc\, the unbinding 
transition of the membrane is continuous. For \U\ > \Utc\, the unbinding tran- 
sition is discontinuous and is then coupled to the phase separation within the 
membrane. At large absolute values \U\ of the sticker binding energy, the sticker 
concentrations of the two coexisting phases vary only slightly with U since the 
majority of the stickers is already bound. 

The important point is that the phase separation occurs at interaction 
strengths \ W\ which are significantly below the critical interaction strength \ Wc\ 
of the 2d lattice gas. At Zy = 0.1, the phase separation can be observed already 
at a fraction W/Wc — 0.31 of the critical mean-field value Wc/kBT — —1 for the 
2d lattice gas. This means that the phase separation is predominantly caused 
by fluctuation-induced interactions between the stickers. 

Monte Carlo simulations confirm this result. Fig. |8] shows a comparison 
of phase diagrams from mean-field theory (lines) and Monte Carlo simulations 
(data points). The cis-interaction strength of the stickers here is W = 0.5Wc. 
In the Monte Carlo simulations, the sticker concentration X = {rii) /a? is deter- 
mined as a function of the chemical potential /i, for various values of the binding 
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energy U . A first-order transition is reflected in a discontinuity of X{^) at a 
certain value of the chemical potential. The two limiting values of X at the 
discontinuity are the concentrations of the coexisting phases. 

The Monte Carlo snapshots of Fig. [9] illustrate the lateral phase separation. 
The snapshots are from canonical Monte Carlo simulations with fixed sticker 
concentration X — 0.45/a^, 0.3/a^, and 0.1/a^. The rescaled potential range 
and cis-interaction strength of the stickers have the same values as in Fig. [8] and 
the sticker binding energy is U = —bksT. The three snapshots thus are taken 
at points within the coexistence region of Fig. |8] Since the overall sticker con- 
centration is kept constant in these simulations, the membrane phase-separates 
into bound domains with high sticker concentration and unbound domains with 
small sticker concentration. The extent of the bound phase shrinks with de- 
creasing sticker concentration. 

The phase behavior strongly depends on the rescaled potential range z^. 
At the short potential range ~ 0.1, phase separation can be observed for 
W/Wc > 0.31 in the mean-field theory, and for W/Wc > 0.35 ± 0.05 in Monte 
Carlo simulations [5]. At the larger potential range z„ = 0.5, phase separation 
occurs for W/Wc > 0.80 in mean-field theory and for W/Wc > 0.65 ± 0.05 
in the simulations. The fiuctuation-induced interactions between stickers thus 
decrease with increasing potential range. The reason for this decrease is that 
bound stickers with larger potential range are less restrictive for membrane 
fluctuations than bound stickers with shorter potential range. 



4.4 Large stickers 

To capture all possible membrane shape fluctuations, the linear patch size of 
the discretized membrane model has to be equal to the cut-off length for the 



fluctuations, see section 2.2 So far, we have only considered stickers with a 
lateral extension which is smaller than or equal to a single membrane patch. 
But the lateral extension of large stickers may be bigger than the cut-off length 
for the membrane ffuctuations. These stickers then occupy several membrane 
patches. 

A simple example are quadratic stickers which occupy Q — 2x2 or 3x3 
lattice sites, see Fig. 10. If each of the Q membrane patches occupied by a sticker 
molecule interacts with the second membrane via a square-well interaction with 
binding energy U and potential range ly, these large stickers can be seen as 
quadratic arrays, or clusters, of Q stickers with the size of a single membrane 
patch. The Hamiltonian of a membrane with these stickers can be written in 
the form [2] 

n{l, n} ^Hei+Y^n,- {V,{1) ~^i)+Y. W^ju^n, (60) 
with the sticker adhesion potential 

Q 

V,{1) = V{k,l, . . . , k,Q) = UY,&ilv - k,q) (61) 

9=1 
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Figure 9: Typical Monte Carlo configurations of a two-component membrane 
containing small, flexible stickers with cis-interaction strength W — O.bWc = 
—ksT ln{l + \/2), binding energy U ~ — Sfc^T and rescaled potential range 
Zy = 0.1. Bound sticker are black, unbound stickers gray, and membrane patches 
without stickers are white. The fixed sticker concentration X is 0.45/a^, 0.3/a^, 
and 0.1/a^ (bottom to top). The membrane phase separates into a bound 
phase with high sticker concentration and an unbound phase with low sticker 
concentration. 
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Figure 10: Membrane segments with: (left) small stickers of size Q = 1 shown 
in black; (middle) quadratic stickers occupying Q — 2 x 2 membrane patches; 
and (right) quadratic stickers occupying (5 = 3x3 membrane patches. 



Here, {(i, 1), . . . , (i, Q)} denotes quadratic arrays of (5 = 2x2 or 3x3 lattice 
sites. The site i = (i,l), for example the central site of a sticker with size 
(5 = 3x3, indicates the sticker position. The cis-interactions in the Hamiltonian 



( 60 ) then are the repulsive hard-square interactions 

Wij = oo for j in Af and zero otherwise (62) 

where Af denotes the 'exclusion area' of an individual (5-sticker at lattice site 
i. The hard-square interactions prevent an overlap of the stickers. 



Phase diagrams for large stickers are shown in Fig. 11 The unbinding tran- 
sition of the membranes is discontinuous for large values of |t/| and continuous 
for \U\ < 2kBT. At large values of \U\, the membrane phase-separates into 
an unbound sticker-poor and a bound sticker-rich phase. The extent of the 
two-phase region increases with the sticker size Q. Since the cis-interactions 



(62) of the stickers are purely repulsive, the phase separation is driven only 
by fluctuation-induced interactions between the stickers. These interactions in- 
crease with the sticker size (see Fig. 11), but decrease with the sticker potential 
range (see Fig. 12). The entropic interactions decrease with increasing potential 
range since bound stickers with larger potential range are less restrictive for the 
membrane shape fluctuations. 

To understand the increase of the fluctuation-induced interactions with the 
sticker size, let us consider an arbritrary shape of the adhering membrane, and 
let us divide the membrane surface into two types of domains: (i) 'Bound do- 
mains' with a membrane separation smaller than the sticker potential range ?„; 
and (ii) 'unbound domains' with a separation larger than 1^. Obviously, in or- 
der to gain more adhesive energy one has to place more stickers into the bound 
domains. If these stickers have the size (5 = 1 of a single membrane patch, the 
maximal adhesive energy depends only on the total area of the bound domains, 
but not on the number of bound domains, or the shape of these domains. In 
contrast, if the stickers are larger and occupy several patches, the total adhesive 
energy depends on the detailed geometry of the bound domains. 

As an example, consider stickers with size (5 = 2x2 and two bound domains 
which both have an area of 2 x 3 patches. As long as these two domains are 
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Figure 11: Phase diagrams for large quadratic stickers occupying (a) Q — 2 x 2 
and (b) Q = 3 X 3 membrane patches. Each of the Q membrane patches 
occupied by a single sticker interacts with the second membrane via a square- 
well potential with binding energy U and rescaled potential range = 0.1. The 
extent of the 2-phase region increases with the sticker size, which indicates an 
increase of the fluctuation-induced interactions between the stickers. 



disjoint, we need four 2x2 sticker to obtain the maximal adhesive energy from 
these two domains. In contrast, if we combine the two bound domains into a 
single 2x6 bound domain, we need only three such stickers in order to gain the 
same adhesive energy. In general, if we cover the area of the bound domains 
with stickers occupying Q > 1 patches, many of these stickers will sit on the 
domain boundaries and thus will not contribute to the adhesive energy in the 
same way as those in the interior of the bound domains. This leads to an 
additional effective line tension of the boundaries between bound an unbound 
domains. This additional line tension favors the aggregation of bound domains 
and corresponds to an increase in the fluctuation-induced interactions. 



4.5 Rigid stickers 

4.5.1 Tensionless membranes 

Stickers with rather stiff 'anchors' in the membrane may change the local mem- 
brane elasticity. If the stickers only affect the local bending rigidity k and do 
not change the modulus of Gaussian curvature, the discretized Hamiltonian has 
the form 3 

i i ^ ^ 

in the absence of cis-interactions between the stickers. Here, Kg is the bending 
rigidity of membrane patches which contain stickers, and k is the rigidity of 
patches without stickers, i.e. the rigidity of the 'bare' lipid bilayer. 
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Figure 12: Phase diagrams for large quadratic stickers with size (5 = 2x2 and 
Q = 3 X 3 as a function of the sticker concentration X and the rescaled potential 
range z„. The binding energy of a single sticker patch has the relatively large 
value U = —bhsT. The sticker concentration in the unbound phase (left data 
points) then is close to zero, see also Fig. 11 The sticker concentration in the 
bound phase (right data points) decreases with increasing sticker potential range 
Zy since the fluctuation-induced interaction becomes weaker. 



Rigid stickers aggregate also without any attractive cis-interactions. The 
fluctuation-induced interactions between rigid stickers thus are signiflcantly en- 
larged. Similar to the case of large sticker considered in the previous section, 
the increased tendency for lateral phase separation can be explained by an ad- 
ditional effective line tension between bound membrane segments [31. Fig. 13 
shows how the phase behavior of the membranes depends on the rigidity Kg 
of the sticker patches. At the rescaled potential range — 0.1 of the sticker 
square-well potential, lateral phase separation occurs for Kg/n > 1.3 ± 0.3 ac- 
cording to Monte Carlo simulations. At Zy — 0.5, the membranes phase-separate 
for Ks/k > 1.5±0.5. This means that the stickers aggregate already at rigidities 
Ks which are only slightly larger than bare lipid bilayer rigidity k. As in the 
previous sections, the fluctuation-induced interactions decrease with increasing 
sticker potential range Zy, see Fig.[l4) 
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Figure 13: Phase diagrams for rigid stickers as a function of the sticker concen- 
tration X and the increased bending rigidity Kg of the membrane patches with 
stickers. The sticker rigidity Kg is given in units of the bending rigidity k for 
the lipid bilayer. The stickers have the binding energy U = —bksT and the 
rescaled potential ranges — 0.1 (left) and z„ = 0.5 (right). The data points 
represent the sticker concentrations in the two coexisting phases, a sticker-poor 
unbound and a sticker-rich bound phase. At the relatively large sticker binding 
energy U — —bhsT considered here, the sticker concentration in the unbound 
phase is very low. 



4.5.2 Effect of tension 

Biological and biomimetic membranes are often under lateral tension. The ten- 
sion suppresses membrane fluctuations and will therefore also affect fluctuation- 
induced interactions. A lateral tension a leads to the additional term 

i 

in the Hamiltonian. Here, 

describes the local area increase of the curved membrane compared to a planar 
membrane configuration with constant separation /. 

Fig. [15] shows how the phase behavior of membranes with rigid stickers de- 
pends on the rescaled lateral tension aa'^/n. The sticker ridigity Kg and binding 



energy U are the same as in Fig. 14 At small lateral tensions a, the concen- 
trations of the coexisting phases agree with the tensionless case, see Fig. [15] 
At higher tensions, the width of the coexistence region decreases since shape 
fluctuations and, thus, fluctuation-induced interactions are suppressed. 

However, the tensions at which the tendency for lateral phase separation 
clearly decreases are relatively high. Lipid membranes typically rupture at ten- 
sions CTr of a few millinewton per meter [6T. Taking cr,. = AmN/m, a membrane 
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Figure 14: Phase diagrams for rigid stickers as a function of the sticker con- 
centration X and the rescaled potential range Zy. The sticker patches have the 
bending rigidity Kg = 1000k, and the sticker binding energy is U = —bhsT as 
in Fig [13] The sticker concentration of the bound phase strongly decreases with 
increasing potential Zy, which indicates a decrease of the fluctuation-induced 
interactions between the stickers. 



patch extension a = 5 nm, and the bending rigidity k = 10~ J leads to the es- 
timate OyC? I K = 1 or logj^Q^arO^ / k) — for the maximum value of the reduced 
tension. The tensions that cause a significant decrease of the coexistence region 
thus are already close to or even above this estimate for the maximum value, 
see Fig. [15] 

To understand this behavior, one has to realize that membrane fluctuations 
are suppressed only on length scales larger than the crossover length ^ nja. On 
smaller scales, thermal fluctuations are still governed by the bending energy. 
The decrease of the coexistence regions in the phase diagrams of Fig. [15] sets 
in at values of the reduced tension aa? /k that correspond to crossover lengths 
of only a few lattice constants. The relevant fluctuations thus turn out to be 
fluctuations of the non-adhesive membrane segments between the small clusters 
of bound stickers, see Fig. [6] A decrease of fluctuation-induced interactions 
then occurs if the crossover length \/ nja is comparable to or smaller than the 
mean distance between the sticker clusters. This interpretation agrees with the 
observation that the influence of an increasing lateral tension is most pronounced 
for the rescaled potential range Zy = 0.5. For this value of 2„, the sticker 
concentration of the bound phase at low tensions is smaller than in the other 
two cases Zy =0.1 and Zy = 0.3, see Fig. [15] The average distance between the 
stickers in the bound phase therefore is larger than in the other cases. 
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Figure 15: Phase diagrams for rigid stickers as a function of the sticker concen- 
tration X and the dimensionless tension ac? jn at different rescaled potential 
ranges z^- As in Fig. [14] the binding energy of the stickers is J7 = —5kBT, 
and the increased bending rigidity of sticker patches is Kg = 1000k. The ex- 
tent of the coexistence regions decreases at large values of a where membrane 
fluctuations on the scale of the average distance between bound stickers are sup- 
pressed. These small-scale fluctuations cause the entropic interactions between 
the stickers. 
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5 Barrier mechanisms for domain formation 



5.1 Membranes with stickers and mobile repellers 

Biological membranes often contain repulsive glycoproteins. These 'repellers' 
form a protective barrier, the glycocalyx. If a membrane contains both stickers 
and repellers, the positions of the molecules can be described by a discrete 
concentration field n with three different values, e.g. the values rii — \ for 
membrane patches i that contain stickers, = 2 for membrane patches with 
repellers, and rii = for patches without stickers or repellers. The grand- 
canonical Hamiltonian of the membrane then can be written in the form [3l [5] 

n{l, n} = Uelil} + J2 I'^l.". (^^(^^) - + {y^k) - Mr)] (66) 

i 

with the Kronecker symbol (5^^ = 1 for i = j and 6ij — otherwise. Here, 
Vs{li) and Vr{li) are the interaction potentials of stickers and repellers, and /is 
and fj,r the chemical potentials for the stickers and repellers. 

Since we have neglected cis-interactions between stickers or repellers, the 
Hamiltonian (66 1 is linear in the concentration field n. Therefore, the degrees 
of freedom of the concentration field can be summed out exactly, see section 
|3.3| This leads to the partition function 



Z = 



[■H.i{l}+T.,V^f{h))]/kE 



of two homogeneous membranes interacting via the effective membrane potential 



1 + Q[t^.-VAh)]/kBT g[M.-K(i,)]/fci 



T 



Here, N denotes the total number of membrane patches. 

In the following, the stickers are characterized by a square-well potential 
with binding energy Us < and range ly, and the repellers by a square-barrier 
potential with barrier energy Ur > and range Ij. > 1^. The range, or 'size', of 
the repellers thus is larger than the sticker binding range. The effective potential 
then has the form 

Vefik) = Uco for 0<li<h 

= Uba for Is <k < Ir 

^ for Ir < k (68) 
with the effective contact energy 

Uco = — fc_Brin Y, — ^ n — ^ (69) 

and the effective barrier energy 

Uba = — fesTln J-. — ^ J-. — — (70) 
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Figure 16: Phase diagrams for a membrane with stickers and repellers as a 
function of the sticker concentration Xs and the rescaled potential ranges and 
Zr of stickers and repellers. The chemical potential for the repellers is /i^ — 0, 
which corresponds to a repeller concentration a^X^ — 0.5 in the unbound phase. 
At large values of z^ and Zr, the repeller barrier with rescaled width 
strong, which leads to lateral phase separation and 2-phase coexistence. The 
data points correspond to the sticker concentrations in the bound phase. The 
sticker contration a^Xg — 10~^ in the unbound phase is given by the vertical 
lines. The critical points obtained from extrapolation are represented by stars. 



The effective barrier energy Uta is positive because of [/,■ > 0, and the effective 
contact energy Uco is smaller than Uba because of Us < 0. Summing out the 
degrees of freedom of the concentration field thus leads to an effective potential 
with a potential minimum Uco induced by the stickers and a potential barrier 
Uha caused by the repellers. 

Bound states of the membrane are only possible if the effective contact energy 
Uco is negative and small enough to compensate the loss of configurational 
entropy, i.e. the entropy difference between the bound and unbound state. This 
entropy difference arises since membrane fluctuations in the bound state are 
more restricted than fluctuations in the unbound state. Suppose the barrier 
energy Uba is high and the membrane is confined to the potential well with 
width Is of the effective potential ( 68 1 . The entropy loss of the bound membrane 
may be estimated via Cfi{kBT)'^ / kI'^ as follows from (56 1 if one replaces the 
average membrane separation I by the width Is of the potential well [68] ■ The 
free energy difference per unit area between the bound and the unbound state 
of the membrane with adhesion potential ( 68 1 can then be estimated as AF = 
— \Uco\+Cfi{kBT)'^ / kI'^ 48J. According to this estimate, the membrane is bound 
for 

\Uco\ > cji{kBTf/nll (71) 

since the free energy difference AF has to be negative to have a stable bound 
state. 

The character of the unbinding transition depends on the strength of the po- 
tential barrier. The barrier induces a line tension between bound and unbound 
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membrane domains, simply because the membrane segments in the boundary 
region between these domains cross the potential barrier and thus have an un- 
favorable potential energy Uta- According to scaling arguments |48l |49], the 
unbinding transition is discontinuous for strong barriers with 

Uba{lr-lsY > ca'ikBTY/K (72) 

and continuous for weak barriers with Ubai^r ~ h)^ < ca'^{kBT)'^ / k. A dis- 
continuous transition implies the coexistence of a bound phase with a high 
concentration of stickers and an unbound phase with a low sticker concentra- 
tion. Sufhciently strong barriers therefore also lead to lateral phase separation 
and sticker aggregation. This barrier mechanism for lateral phase separation is 
weaker at higher temperatures T, in contrast to the entropic mechanisms dis- 
cussed in chapter 4. Higher temperatures T require larger barriers Uba for phase 



separation, see eq. (72 1 



The coefficient c in eq. ( 72 1 can be estimated from Monte Carlo simulations 
[5]. In the simulations, the sticker concentration Xg = i)/a^ is determined 
as a function of the chemical potential fig of the stickers. Lateral phase separa- 
tion is reflected in a discontinuity of Xs{fis) at a transition value fig = fJ-t- The 
two limiting values of Xs at fi* correspond to the sticker concentrations in the 
two coexisting phases. Monte Carlo phase diagrams obtained from such simu- 



lations are shown in Fig. 16 The coefficient c in eq. (72) can be estimated from 
the critical points for the 2- phase coexistence regions of the diagrams. These 
critical points are at = 0.15 ± 0.05 in the left and at 2^ = 0.65 ± 0.05 in 



the right diagram. In agreement with eq. (72), the two values differ by a factor 
of 5 within the numerical accuracy, and thus reflect the same critical barrier 
strength with c = 0.013 ± 0.005. 



5.2 Membranes with stickers and generic repulsive inter- 
actions 

In the last section, we have considered the interplay of stickers and repellers 
during membrane adhesion. A related situation arises if stickers act against a 
generic repulsive interaction potential of the membranes. 'Generic' means that 
the repulsive interaction is, or is taken to be, independent of the local compo- 
sition of the membrane, in contrast to the 'specific' sticker interactions. One 
example for such a generic interaction is the electrostatic repulsion of equally 
charged membranes, provided the charge distributions can be approximated as 
uniform distributions. Another example is a repulsive brush of immobilized re- 
peller molecules. For cell membranes, immobilization of membrane molecules 
can arise from an anchoring to the cytoskeleton of the cell. Immobile repellers 
have approximately fixed positions in the membrane, in contrast to the mobile 
repellers considered in the last section which are free to diffuse within the plane 
of the membrane. 

A discretized membrane with stickers and a generic potential Vg{li) can be 
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described by the Hamiltonian [1] 



i 

if cis-interactions between the stickers are neghgible. Summing out the sticker 
degrees of freedom in the partition function Z as described in section 3.3 then 
leads to 



z = 



n 



exp 



with the effective potential 

Vef{l)^Vg{l)-kBT\n\l+c^^ 



- Vs{l) 
knT 



(74) 



(75) 



5.2.1 Stickers with square-well potential 

Let us first consider stickers which are again characterized by a square-well 
potential Vs{li) = U6{ls — li) with binding energy U and range /g, and a repulsive 
generic interaction which is characterized by a square-barrier potential Vg{li) = 
Ubadiha — h) with barrier energy Uba > and range Ir- If the range Ir of the 



generic repulsion is larger than the sticker range 1^ , the effective potential ( 75 1 
has the form 

Vefik) - K = Uco for 0<k<ls 

= Uba for Is < k < Ir 

= for Ir < k (76) 

with the contact energy 

lj^g{l^-U)/kBT 

Uco = Uba - keT In ^ ^ ^^^^^^ (77) 

The constant term Vo — — fesTln (l + e^/*^^-^) depends only on the reduced 
chemical potential fi/ksT of the stickers. 



This potential has the same form as the effective potential ( 68 1 for a mem- 
brane with stickers and (mobile) repellers, which implies also the same phase 
behavior: The membrane is bound if the contact energy \Uco \ exceeds the thresh- 



old value estimated in eq. ( 71 1 , and phase-separates if the barrier exceeds the 



critical strength estimated in eq. ( 72 ) 



5.2.2 Stickers with linear potential 

In the previous sections, we have seen that a barrier in the effective membrane 
potential causes lateral phase separation. This barrier mechanism for phase sep- 
aration is rather general. So far, we have considered the interplay of attractive 
square-well and repulsive square-barrier potentials. For a deeper understanding 
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of the mechanism, it is instructive to consider also the analytically convenient 
case in which the generic potential can be approximated by a harmonic potential 



V2_ 

2a2 



(78) 



If a generic potential Vg has a relatively deep minimum at a certain separation 
lo of the membranes, the harmonic approximation (78 1 can be justified by a 
Taylor expansion around the minimum. The prefactor in eq. ( 78 ) is then V2 = 
a\<fVjdP)\i^^. 

Let us further assume that the corresponding sticker potential Vs{l) has 
an essentially constant gradient for those values of I for which we can use the 
harmonic approximation ( 78 ) for the generic potential. In such a situation, we 
may truncate the expansion of the sticker potential in powers of I ~ and use 

ma 

a{l - lo) 



V,{1) = VSo) 



(79) 



with a = adVs{l)/dl\i^ > 0. This approximation may be valid for extensible, 
or 'spring-like', sticker molecules that are irreversibly bound to the membranes 
and have an unstretched extension small compared to lo- 

To simplify the notation, we introduce here the dimensionless parameters 



a 



and the rescaled separation field 

Zi = a/2 + 



and 



M - 



Vs{lo) 



V2 h - lo 



(80) 



(81) 



Let us first consider the special line in the (/i, a) parameter space given by 

/2 = /i* = -aV2 , (82) 



Along this line, the effective potential resulting from the eqs. (751, (78 1 and ( |79| 
can be written in the form 



Vef{z) 

ksT 



-ln[2cosh(az/2)] 



(83) 



which is symmetric under the inversion z ^ —z. When one varies the param- 
eter a while keeping /i = /i*(a), the effective potential exhibits a continuous 
bifurcation at the critical value a = etc = 2, see Fig. 17 The potential has a 
single minimum for a < etc, and two degenerate minima for a > cxc- The critical 
value CKc — 2 of the bifurcation point can be directly inferred from the second 
derivative of eq. ( 83 1 : 



1 d^V.fiz) 



ksT dz^ 



4cosh^(az/2) 



(84) 
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Figure 17: The effective potential Vef as a function of tlie rescaled separation z 
for three values of the coupling parameter a. The analytical expression for Vef 
is given in (83 1. The effective potential has a single minimum at a = 1.5 (top), 
and two degenerate minima at a — 2.5 (bottom). At a — etc — (middle), the 
potential undergoes a continuous bifurcation. 
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For z — 0, this expression is equal to 1 — (5^/4, which vanishes for a = fie = 2. 

Limit of rigid membranes - At large values of the ratio k/v2, the membrane 
is quasi rigid. Thermally excited shape fluctuations of the membrane thus can 
be neglected. The free energy T = —{kBT/A)\uZ per area A is then simply 
given by Vg^f/a?, and the phase behavior can be determined by minimizing the 



effective potential. For jl = il»{a) and a > 2, the effective potential (83 1 is a 
symmetric double-well potential with two degenerate minima. As soon as the 
chemical potential jl deviates from its coexistence value = jla., this symmetry 
is broken and the effective potential has a unique global minimum. The system 
thus exhibits a discontinuous transition when one changes the chemical potential 
from /i = /i* — e to /i = /i* + e for a > 2. Here, e denotes a small chemical 
potential difference. In the limit of rigid membranes, the critical point for phase 
separation is identical with the bifurcation point of the effective potential at 
etc = 2 and jlc = —<^c 12 — —2. 

The positions of the extrema of the effective potential are determined by 
dVg^fi^z) I dz = 0. Along the coexistence line with /i = /i* = —a^/2, this leads to 
the transcendental equation 

a , / cy.z\ 
z--tanh(^ — j . (85) 

This equation has the trivial solution z — Q for all values of a. This solution 
corresponds to a minimum for a < olc — 2, and to a maximum for a > oic = 2. 



For a > = 2, eq. (85) has two additional solutions corresponding to the two 



degenerate minima of the effective potential 14/, see Fig. 17 
Along the coexistence line, the sticker concentration 



is given by 

X 



(87) 



The concentrations of the coexisting phases then are obtained by inserting the 



numerical solutions of the transcendental equation (85 1 into (87). The resulting 
phase diagram is shown in Fig. |18[ Inside the shaded two-phase region, a sticker- 
poor phase with large membrane separation coexists with a sticker-rich phase 
with smaller separation. 

Flexible membranes - A flexible and, thus, fluctuating membrane can easily 
cross small barriers in the potential. First-order transitions then only occur 
if the barrier exceeds a critical height |351 For a flexible membrane, the 
critical coupling constant etc therefore will be larger than the bifurcation value 



cic = 2 of the effective potential ( 83 1 



With Monte Carlo simulations, the critical point can be determined via the 
moments ^ ^ 

C2 — ,|_|,, and C4 — ,_„, „ (88) 
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Figure 18: Phase diagram for linear stickers in the absence of membrane fluctu- 
ations, depending on the sticker concentration X and the reduced couphng con- 
stant a. Within the grey two-phase region, a sticker-poor phase characterized 
by a relatively large membrane-surface separation coexists with a sticker-rich 
phase for which this separation is relatively small. The critical point is located 
at a'^Xc — 1/2 and = 2. 



Here, 



1 ^ 

i=l 



(89) 



is the spatially averaged order parameter, and (• • • ) denotes averages over all 
membrane configurations |49l 154] . In principle, the values of these moments de- 
pend on the correlation length ^ and the linear size L of a membrane segment. 
But at the critical point, the correlation length ^ diverges, and the values of the 
moments become independent of L [151 ISl]. Therefore, the critical coupling con- 
stant fie can be estimated from the common intersection points of the functions 



C2(fi) and C4(fi) at different values of L, see Fig. 19 



Fig. [20] shows the critical rescaled coupling constant fic as a function of the 
reduced rigidity k/v2. For large k/v2, etc approaches the limiting value fic = 2 
of 'rigid' membranes, see above. With decreasing k/v2, the membrane shape 
fluctuations become more pronounced and lead to an increase in the value of 
fic- As in section 5.1 the membrane fluctuations thus reduce the tendency for 
lateral phase separation. 

Lateral phase separation occurs for coupling constants fi > fic. In either of 
the two phases, the entire membrane then is located around one of the minima 
in the effective potential. In the sticker-poor phase, the membrane is found in 
the minimum with larger membrane separation. This minimum is dominated 
by the generic membrane potential and corresponds to a state of weak adhesion. 
In the sticker-rich phase, the membrane fluctuates around the minimum with 
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Figure 19: Monte Carlo data for the moments C2 and C4 defined in eq. (88 1 as 
a function of the reduced coupHng constant a. The ratio of the bending rigidity 
K and the strength vi of the generic harmonic potential ( 78 1 has the fixed value 
k/v2 = 1. The membrane segments considered in the simulations consist of 
N X N membrane patches. 



smaller separation, which corresponds to a state of tight adhesion. In contrast, 
there is only a single phase for coupling constants a < ac- For 2 < a < etc for 
example, the two minima of the effective potential are both populated by many 
different segments of the fluctuating membrane, which is then able to cross the 
potential barrier between the minima. 

6 Dynamics of domain formation during adhe- 
sion 

So far, we have focused on equilibrium aspects of the domain formation. In 
this section, we consider the adhesion dynamics, or in other words, the time- 
dependent evolution of the domain patterns. The models presented in this 
section mimick the adhesion geometry of cells or vesicles by dividing the mcm- 
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Figure 20: The critical values Sc of the coupling constant a as a function of 
k/v2- Here, k is the membrane rigidity, and V2 is the strength of the generic 
harmonic potential. The coupling constant a is defined in eq. ( 80 ) and governs 



the strength of the linear sticker potential (79 1. For large values of k/w2, the 
critical coupling constants tend towards the value cfc = 2 of 'rigid' membranes. 
The statistical errors here are smaller than the symbol sizes. 



branes into two zones: a contact zone, and a surrounding nonadhering mem- 
brane region. The stickers and receptors can diffuse in the whole membrane, but 
interact with the second membrane only within the contact zone. The problem 
of modeling the full shapes of cells or vesicles thus is avoided in these models. 
Instead, the contact zone is assumed to have an essentially circular shape and 



a constant area on the relevant time scales, see Fig. 21 For biomimetic vesicles 
with stickers and repellers, this adhesion geometry can result from fast initial 
gravity-induced adhesion of the vesicles on a supported membrane |69j . In the 
case of adhering T cells, initial adhesion is mediated by relatively long integrins. 
Experimental pictures show that the contact zone of T cells fully develops in 
less than 30 seconds [501. 



6.1 Adhesion of vesicles with stickers and repellers 

Within the contact zone, the membrane of an adhering vesicle with stickers and 
repellers can be described by the Hamiltonian [5j 



(90) 



where Si,j is the Kronecker delta. In the nonadhering membrane region sur- 



rounding the contact zone (see Fig. 21 1, the stickers and repellers do not in- 
teract with the second membrane. In this region, the configurational energy of 
the stickers and repellers thus is constant, i.e. independent of the sticker and re- 
peller positions, and independent of the membrane shape. The membrane shape 
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Figure 21: 'Cell' adhesion geometry: The circular contact zone is surrounded 
by a nonadhering membrane ring. Stickers, or receptors, diffuse around in the 
whole membrane, but interact with the apposing membrane only within the 
contact zone. 



therefore is only modeled explicitly within the contact zone via the separation 
field /. In this model, we use 'free' boundary conditions of the separation field / 
at the contact zone rim, i.e. an unconstrained boundary separation. The Hamil- 
tonian ( 90 1 is the canonical equivalent of the grand canonical Hamiltonian ( 66 1 



of section 15.11 We consider now the canonical ensemble since the total numbers 
of stickers and repellers in the whole vesicle membrane are constant. 

In the following, the sticker potential is a square-well potential with depth 
Us = —lOkgT and rescaled range Zg, and the repeller potential is a square- 
barrier with height Ur = lOfc^T and rescaled range Zr > Zg. For these rela- 
tively large energies of lOfc^T, the sticker binding is nearly irreversible, i.e. the 
majority of stickers stays bound after first binding, and the repellers exclude 
membrane separations smaller than almost completely. The contact zone of 
the membrane here has the diameter 100 a, and the whole membrane has the 
diameter 200 a, see Fig. [21] 

The formation and evolution of the domain patterns can be studied with 
Monte Carlo simulations. The Monte Carlo simulations presented below start 
from a random distribution of stickers and repellers, and a rescaled membrane 
separation Zi = Zr in the contact zone, with all stickers unbound. A Monte 
Carlo step consists in attempts (i) to move each sticker and repeller to one of 
the 8 neighbor sites (lateral diffusion) and (ii) to shift the rescaled membrane 
separation Zi at every lattice site i in the contact zone (shape fluctuations). 

In this model, three different dynamic regimes of pattern formation can be 
observed. The dynamic regimes depend on the characteristic lengths and the 
concentrations of stickers and repellers: 

Regime (A): Long repellers impose a strong barrier to sticker adhesion. The 
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Figure 22: Typical pattern evolution in the dynamic regime (B). The potential 
ranges of stickers and repellers are Zg — 1.0, Zr — 1.6 and the overall concentra- 
tions of stickers and repellers are a?Xs = 0.1/a^, a^X^ = 0.5/a^. Stickers are 
shown in black, repellers in grey. The circle represents the rim of the contact 
zone. Because of the diffusion of stickers into the contact zone, clusters at the 
rim grow faster, forming rings at intermediate time scales. The final configu- 
ration represents the equilibrium state. The snaphots are taken at 10^, 4 • 10^, 
lO'', 4 • 10^, 1.6 • 10^ 6.3 • 10^ 10^ and 4 • 10^ MC steps. 



nucleation time for sticker binding therefore is large compared to typical diffu- 
sion times, and the membrane binds via growth of a single sticker nucleus. Such 
an adhesion behavior has been recently observed for biomimetic vesicles with 
PEG-lipopolymers as repellers and integrins as stickers (69j. 

Regime (B): For short repellers, the nucleation time for sticker binding is 
small, and many nuclei of bound stickers are formed initially. Since unbound 
stickers diffuse into the contact zone, nuclei at the rim of this zone grow faster. 
At sufficiently high sticker concentrations, an intermediate ring of bound stick- 
ers is then formed, enclosing a central domain of repellers, see Fig. |22[ Later, 
this pattern inverts, and a single central sticker cluster is surrounded by re- 
pellers. The central position of this sticker clusters here is caused by membrane 
shape fluctuation in the unbound repeller domain surrounding the cluster. The 
sequence of patterns in this regime has a striking similarity to the pattern evo- 
lution observed during T cell adhesion (see Introduction and next section). 

Regime (C): In an intermediate regime, the sticker concentration is not large 
enough for the formation of a closed sticker ring from the initial nuclei. Instead, 
circular arrangements of separate sticker clusters emerge at intermediate times, 
see Fig. [23] 

Two quantities are helpful to characterize the three different regimes (A), 
(B), and (C) systematically [5]. The first quantity is the maximal number 
of bound sticker clusters _/V^°^ during the pattern evolution. The number of 
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Figure 23: Typical pattern evolution in the dynamic regime (C). The potential 
ranges of stickers and repellers are Zg — 1.0, Zr = 1.8, and the overall concen- 
trations are Xg = 0.04/a^, = 0.5/a^. Stickers are shown in black, repellers 
in grey. The sticker concentration is smaller than in Fig. |22] Therefore, circular 
patterns of separate clusters emerge at intermediate times, instead of closed 



sticker rings as in Fig. 22 The snaphots are taken at 10^, 4 • 10^, 1.6 • 10^ 



6.3 • 10*, 4 • 10^ 6.3 • 10^ 1.6 • 10^ and 4 • 10^ MC steps. 



sticker clusters Nd first increases since new sticker clusters nucleate, and later 
decreases as a consequence of cluster coalescence. The second quantity is the 
maximal sticker occupation in an outer ring of the contact zone. To 

define this quantity, the membrane ring with distances 40 a < r < 50 a from 
the center of the contact zone is divided into 100 equal segments, each covering 
an angle 27r/100. The ring occupations simply is the fraction of segments Y 
which contain bound stickers. The ring occupation Y has a maximal value 
ymax intermediate times when a ring of sticker clusters is formed. 

Appropriate values to describe the crossover between the three dynamic 
regimes are the maximal number of sticker clusters iV™"^ = 3 and the maximal 
ring occupation Y"^"-^ = 0.8 (see Fig. 24): Simulations with Y™"-^ > 0.8 show 
intermediate configurations with a closed ring of bound stickers as in Fig. |22| 
(dynamic regime (B)). For iV™°^ < 3, on the other hand, adhesion proceeds by 
sticker condensation mostly around a single, dominant nucleus. For 7V™"^ > 3 
and y""^^ < 0.8, configurations with a circular arrangement of separate clusters 
emerge as in Fig. 23 (dynamic regime (C)). 
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Figure 24: Dependence of the dynamic regimes (A), (B), and (C) on the sticker 
and repeller concentrations Xg and Xj. and the repeller range z^. The sticker 
potential range here is Zg = 1.0. The repeller concentration is Xj. = 0.5/a^ in the 
left diagram. In the right diagram, the repeller range is z^ — 2.0. Characteristic 
intermediate patterns are indicated for all three regimes. In the upper left grey 
area (regime (A)), the maximal number of sticker clusters iV™"^ during adhesion 
is smaller than 3. In the lower right grey area (regime (B)), the maximal ring 
occupation y""*^ is larger than 0.8. The nucleation time for sticker clusters 
increases with the length Zr or concentration Xr of repellers, and thus also the 
extent of the dynamic regime (A). 



6.2 Adhesion of T cells 
6.2.1 Model 

Helper T cells mediate immune responses to antigen-presenting cells ( APCs) , see 
Introduction. The cell adhesion model presented here considers two apposing 
membranes with different concentration fields. The first membrane represents 
the T cell and contains the receptors TCR and LFA-1. The second membrane 



represents the APC and contains the ligands MHCp and ICAM-1 (see Fig. 25 1. 

To keep MC simulations of the relatively large cell membranes and contact 
zones tractable, the value z = 1 of the rescaled membrane separation field 
z — {I / a)yj K./T corresponds in this model to a length of 20 nm. This results 
in the relation a — 20-y/ K/{kBT) nm for the linear patch size. For the typical 
bending rigidities ki — K2 — 25A;bT' of the two biomembranes ^70], the effective 
rigidity k has the value VZ.bksT, and the linear patch size is a ~ 70 nm. Monte 
Carlo simulations with smaller patch sizes should lead to comparable results, 
but require significantly longer computation times. 

Since the membrane patch size a is relatively large, the model allows several 
molecules in a single patch. The local composition of the T cell membrane is 
then described by the numbers nf of TCRs, nf of LFA-1, and nf*" of glyco- 
proteins in each membrane patch i. Correspondingly, the composition of the 
APC membrane is given by the numbers of MHCp, nl of ICAM-1, and 
nf " of glycoproteins in all patches. The overall configurational energy of the 
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Figure 25: Cartoon of a T cell membrane (top) adhering to an APC membrane 
(bottom). The T cell membrane contains the T cell receptor TCR (green) 
and the receptor LFA-1 (red). The APC membrane contains the corresponding 
ligands MHCp (green) and ICAM-1 (red). Both membranes contain repulsive 
glycoproteins (grey). Because of the different lengths of bound TCR/MHCp 
complexes, LFA-l/ICAM-1 complexes, and glycoproteins, the membrane phase 
separates into domains. 



membranes in the contact zone is the sum of the elastic energy ( 92 1 and the 
interaction energies of receptors, ligands, and glycoproteins [B]: 

n{l,7i}=n,i{l} + J2[^HnI,nf)VTM{k)+i-nHnf,ni)Vu{h) 



in?' + n?'') 



Vcih) (91) 



Here, VtmCi) and VLi(^i) are the attractive interaction potentials of TCR/MHCp 
and LFA-l/ICAM-1 complexes, Vcili) is the repulsive interaction potential of 
the glycoproteins. The term min(n^, nf^) denotes the minimum of the numbers 
of TCR and MHCp molecules at site i. This minimum is equivalent to the 
number of interacting TCR/MHCp pairs in the apposing patches at site i. The 
elastic energy 

^ei{?} - [('«/2fl')(Ad^.)' + {cj/2){Vah?] (92) 

i 

of the model has a contribution from a lateral tension a. In the simulations 
presented here, the value for the lateral tension is u = Q.\n/a^ ~ 2 • lO^^A^/m, 
which is within the range of values measured for Dictyostelium discoideum cells 

The receptor-complexes can only form if the membrane separation is in an 
appropriate range. The length of the TCR/MHCp complexes is about 15 nm, 
while the LFA-l/ICAM-1 complexes have a length of about 40 nm [5S]. The 
membrane within a patch is 'rough' because of the thermal fluctuations on length 
scales smaller than the linear extension a ~ 70 nm of the patches. Therefore, re- 
ceptor/ligand complexes can also form if the separation of two apposing patches 
deviates slightly from the precise lengths ztm and zn of the complexes. In the 
model, the interaction potential of TCR and MHCp is characterized by the 



52 



square-well potential 



Vtm 



= 



Utm for 10 nm < k < 20 nm, 
otherwise 



(93) 



and the interaction potential of ICAM-1 and LFA-1 is 



Vli 



U LI for 35 nm < < 45 nm, 
otherwise 



= 



(94) 



Here, Utm > is the binding energy of a TCR/MHCp complex, and Uli > 
the binding energy of LFA-l/ICAM-1. 

The repulsive glycoproteins protruding from both membranes vary in size. 
However, many of these proteins have a length comparable to the length of 
the LFA-1 /ICAM-1 complexes. These glycoproteins do not inhibit the bind- 
ing of ICAM-1 and LFA-1, but impose a steric barrier for the formation of 
TCR/MHCp complexes. They are characterized by the potential 



with Ug ~ IOk/o^ and = 40 nm. This potential aims to capture that a 
membrane patch of size a containing a glycoprotein has to bend around this 
protein to achieve an overall patch separation smaller than the length of the 
glycoprotein. 

6.2.2 Adhesion dynamics without cytoskeletal transport 

We first consider the pattern formation in the absence of active forces which 
transport molecules in or out of the contact zone. In the absence of active 
transport, the diffusive motion of the macromolecules is modeled as an unbiased 
hopping process between neighboring membrane patches as in section |6.1[ Each 
receptor, ligand, or glycoprotein in a certain membrane patch can hop to one of 
the four nearest neighbor patches during a single time step. During a time step, 
we also attempt to shift the separation li between apposing membrane patches 
in the contact zone by d • C[— 1, 1] where d is the step width 10 nm, and C[— 1, 1] 
is a random number between —1 and 1. 

A single Monte Carlo step roughly corresponds to 1 ms of real time. This 
time estimate can be derived from the 2-dimensional diffusion law (x^) — ADt 
and the typical diffusion constant 13 ~ 1 /im^/s for membrane-anchored macro- 
molecules. In a single Monte Carlo step, a free receptor, free ligand, or a glyco- 
protein moves a distance a to a neighboring membrane patch, which corresponds 
to a diffusion time t — /{AD) ~ 1 ms for a = 70 nm. On the length scale 
of our patches, the diffusive motion of the macromolecules is slower than the 
relaxation of the membrane separation [72j and hence defines the time scale. 

As initial conformation of the MC simulations, the separation profile is I = 
lo + cr"^ where r is the distance from the center of the contact zone, Iq is 45 nm. 



Vg 



Ug{1 - for I < Ig, 



otherwise 



(95) 
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Figure 26: Typical pattern evolution without active TCR transport in the dy- 
namic regime 1. Membrane patches with bound TCR/MHCp complexes are 
shown in green, patches with bound LFA-l/ICAM-1 complexes in red. The 
black circle represents the contact zone rim. The effective binding energies 
of the TCR/MHCp and the LFA-l/fCAM-1 complexes are Utm = e.SfcsT 
and Uli = 3A:bT. The overall concentrations of TCR, ICAM-1, LFA-1, and 
glycoproteins in each of the membranes is 0.4/a^ ~ 80 molecules/ /im^ for a 
linear patch size a ~ 70 nm, and the concentration of MHCp is 0.1/a^ ~ 20 
molecules / ^.m^ . 



and c > is chosen so that the separation at the rim of the contact zone with 
radius r — Aba is 100 nm ('clamped' boundary condition). This initial separation 
in the contact zone is larger than 45 nm, and thus beyond the interaction range 
of receptors, ligands, and glycoproteins. Initially, these molecules are taken to 
be randomly distributed within the whole membrane. 

Since the length difference of the complexes leads to phase separation at the 
molecular concentrations considered below, the two types of receptor/ligand 
complexes have to 'compete' for the contact zone. In general, the overall 
area of TCR/MHCp domains in the contact zone increases with the concen- 
trations of TCR and MHCp molecules and with the effective binding energy 
Utm- However, if the molecular concentrations or the binding energy are too 
small, TCR/MHCp domains do not form, and the contact zone contains only 
bound LFA-1 /ICAM-1 complexes. At molecular concentrations and binding en- 
ergies where TCR/MHCp and LFA-l/ICAM-1 domains coexist, we observe two 
different regimes for the dynamics with clearly distinct patterns of TCR/MHCp 
domains at intermediate times. The pattern evolution roughly depends on the 
overall area of TCR/MHCp domains after initial relaxation. 

Regime 1: If the overall area of TCR/MHCp domains is relatively large, 
we observe a characteristic ring-shaped TCR/MHCp domain at intermediate 



times, surrounding a central domain of LFA-l/ICAM-1 complexes, see Fig. 26 



The ring finally breaks to form a single large TCR/MHCp domain. The pattern 



evolution in this regime is very similar to regime (B) of section 6.1 



Regime 2: For smaller TCR or MHCp concentrations, or smaller effective 
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Figure 27: Typical pattern evolution without active TCR transport in the dy- 
namic regime 1. Membrane patches with bound TCR/MHCp complexes are 
shown in green, patches with bound LFA-l/ICAM-1 complexes in red. The ef- 
fective binding energies are Utm — b.bksT and Uli = ^^ksT, and the molecular 
concentrations are the same as in Fig. 



26 



binding energy, characteristic multifocal TCR/MHCp patterns emerge at inter- 
mediate times. A typical example is shown in Fig. |27[ The patterns in this 



regime are similar to those formed in regime (C) of section 6.1 



To distinguish the two dynamic regimes systematically, it is useful to consider 
again the maximal sticker occupation Y"^^^ in a peripheral ring of the contact 
zone, now with distances 35 a < r < 45 a from the center. A pattern evolution 
with Y"^^^ < 80 % typically has multifocal intermediates as in Fig. 27 (Regime 
2) , while pattern evolutions with Y"^'^'^ > 80 exhibit the inverted synapse of T 



cells with a peripheral TCR/MHCp ring, see Fig. 26 (Regime 1) 



The left diagram in Fig. [28| shows how the dynamic regimes for pattern for- 
mation depend on the effective binding energies Utm and Ulj of the TCR/MHCp 
and LFA-l/ICAM-1 complexes. An increase in Utm in general leads to more 
TCR/MHCp complexes in the contact zone, while an increase in Uli leads 
to the binding of more LFA-l/ICAM-1 complexes. Instead of varying the ef- 
fective binding energies Utm and Uli, the numbers of bound receptor/ligand 
complexes in the contact zone could also be changed by varying the overall 
concentrations of the receptors and ligands, with similar effects on the pattern 
formation. 



The right diagram in Fig. 28 shows the effect of the glycoprotein concentra- 
tion Xq on the adhesion dynamics. The length of the glycoproteins is compati- 
ble with the length of the LFA-l/ICAM-1 complexes. Hence, the glycoproteins 
can enter the red LFA-l/ICAM-1 domains in the contact zone, but are excluded 
from the green TCR/MHCp domains. The accessible membrane area for the 
glycoproteins increases with the fraction of LFA-l/ICAM-1 domains in the con- 
tact zone, and so does the entropy of the glycoprotein distribution. Therefore, 
an increase in the overall glycoprotein concentrations leads to a larger fraction 
of red LFA-l/ICAM-1 domains in the contact zone, and thus has a similar effect 
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Figure 28: Dynamic regimes for T cell adhesion. At large values of the bind- 
ing energy Utm of the TCR/MHCp complexes, a peripheral TCR/MHCp ring 
emerges at intermediated times as in Fig. 26 (Regime 1). At medium values 
of Utmj multifocal patterns as in Fig. 27 arise at intermediate times (Regime 
2). At small values of Utm, TCR/MHCp domains in the contact zone do 
not form. The threshold for the formation of TCR/MHCp domains and the 
crossover between the two dynamic regimes depend on the binding energy Uli 
of LFA-l/ICAM-1 complexes and the glycoprotein concentration Xq in both 
membranes. The concentrations of TCR, LFA-1, and ICAM-1 are 0.4/a^ ~ 80 
molecules//xm^ and the concentration of MHCp is 0.1 /a^ ~ 20 molecules//im^ . 
In the left diagram, the glycoprotein concentration in each of the membranes is 
Xg — 0.4/a^. In the right diagram, the binding energy Uli of LFA-l/ICAM-1 
complexes has the value SksT. The black diamonds in the figure represent data 
points obtained from Monte Carlo simulations. 



as increasing the binding energy Uli of the LFA-l/ICAM-1 complexes. 

In both dynamic regimes of pattern formation, the coalescence of clusters 
finally leads to a single TCR/MHCp domain in our model. In the absence of 
active transport processes, we always observe that the final TCR/MHCp domain 
is in contact with the rim of the contact zone, see Figs. [26] and [27] This behavior 
can be understood from the line tensions at the domain boundaries and at the 
rim of the contact zone [6^. In these equilibrium conformations, the length of 
the boundary between the two domains is significantly shorter than in the case 
of a central circular TCR/MHCp cluster. 

6.2.3 Adhesion dynamics with active transport of TCRs 

In T cells, active processes transport receptors into the contact zone [73 and 
glycoproteins out of this region '74' '75] . The framework enabling these transport 
processes is the actin cytoskeleton, which polarizes during adhesion, with a focal 
point in the center of the contact zone [Tj [7B]. For TCRs, the transport is 
mediated by myosin, a molecular motor protein binding to the actin filaments. 
The transport of TCRs can be modeled as a directed diffusion. In the model, 
each TCR molecule is simply assumed to experience a constant force which is 
directed towards the center of the contact zone midpoint. This force corresponds 
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Figure 29: Pattern evolution with active transport of TCRs towards the center 
of the contact zone. Membrane patches with bound TCR/MHCp complexes 
are shown in green, patches with LFA-l/ICAM-1 complexes in red. Molecular 
concentrations and binding energies are the same as in Fig. 26 (dynamic regime 
1). (Top) At zero force, the intermediate TCR/MHCp pattern is stable for 
30 minutes and more. In the final equilibrium pattern, both types of domains 
are in contact with the rim of the adhesion region, see section III. (Middle) 
At the force F = O.OlfcsT/a, the final equilibrium state is the target-shaped 
mature synapse of T cells. This state is already established within 30 minutes. 
(Bottom) At the 10-fold stronger force F = 0.1 ksT/a, the final target-shaped 
pattern already forms within 5 to 10 minutes. An intermediate pattern with a 
TCR/MHCp ring appears around 30 seconds after initial contact. 



to an additional term _F • r in the configurational energy of each TCR where F 
is the magnitude of the force and r the distance of the receptor from the center 
of the contact zone. Under the influence of this force, diffusive steps bringing 
TCRs closer to the focal point of the cytoskeleton in the center of the contact 
are, in general, more likely than diffusive steps in the opposite direction. 

Figure [29] compares the pattern evolution at zero force with patterns at the 
forces F = OmkeT/a ~ 6 • lO^^^ N and = O.lfcsT/a ~ 6 • lO^^^ N. The con- 



centrations and binding energies are the same as in Fig. 26 For these values, 
the force F — O-Olfc^T/a is close to the force threshold leading to a target- 
shaped flnal synapse with central TCR/MHCp cluster. Besides leading to a 
central TCR/MHCp cluster, the active forces speed up the pattern evolution. 
At the weaker force F = O.OlfcsT/a, the flnal equilibrium state is reached af- 
ter approximately 30 minutes, while the 10-fold stronger force F = Q.lksT/a 
leads to equilibrium within few minutes. The absolute times are based on the 
estimate that one Monte Carlo step roughly corresponds to 1 ms, see above. A 
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TCR/MHCp ring at intermediate times is formed in all three cases shown in 
Fig. [29] The forces in the model are average forces acting on a single TCR. 
Since the transport of a TCR molecule over larger distances presumably in- 
volves several cytoskeletal binding and unbinding events, these average forces 
are significantly smaller than the local maximum forces around 1 pN= 10~^^ N 
which can be exerted by a single molecular motor |77j . 

Experimentally, the mature synapse of T cells has been observed to form 
on timescales between 5 and 30 minutes [30l [33] . These timescales agree with 
the equilibration times for the force range O.Ol/csT/a ^ F < OAksT/a with 



T-cell like pattern evolution, see Figs. 29 In the absence of active forces {F 



0), the intermediate peripheral TCR/MHCp ring seems to be metastable and 
appears in our simulations for times up to an hour. This metastability might 
explain the inverted NK cell synapse which consists of a peripheral ring of short 
receptor/ligand complexes, and a central domain containing the longer integrins. 
The inverted synapse of NK cells seems to be formed by self-assembly, since it 
is not affected by ATP depletion or cytoskeletal inhibitors [75] [73] . 

Some groups [BUI ISIl [H21 E3\ have proposed that the final T cell pattern 
can be obtained by self-assembly. In the model of Qi et al. [50], the central 
TCR/MHCp domain apparently results from the circular symmetry of the con- 
sidered patterns. This symmetry prevents patterns with a single TCR/MHCp 
domain at the contact zone rim. Coombs et al. [S3] investigate equilibrium as- 
pects of T cell adhesion and focus on circular symmetric patterns similar to Qi 
et al. In the models of S.-J. Lee et al. [51] and Raychaudhuri et al. [52], the 
central TCR/MHCp domain seems to arise from the boundary condition that 
the membrane separation at the contact zone rim is close to the LFA-l/ICAM-1 
length of 40 nm. This boundary condition favors LFA-l/ICAM-1 domains at 
the rim, and repels TCR/MHCp domains from the contact zone rim. How- 
ever, directly adjacent to the contact zone of two cells, the membrane separa- 
tion quickly attains values much larger than the lengths of the receptor/ligand 
complexes. Therefore, it seems more realistic to impose a membrane separa- 
tion significantly larger than the lengths of LFA-l/ICAM-1 and TCR/MHCp 
complexes as boundary at the rim. In the discrete model presented here, the 
separation at the rim is 100 nm. 

Active cytoskeletal processes may also play an important role for the multi- 
focal patterns of adhering thymozytes [41]. These patterns resemble the char- 
acteristic intermediate domain patterns in the dynamic regime 2, see Fig. [27] 
However, the intermediate patterns are only stable on the timescale of minutes. 
After a few minutes, domain coalescence leads to a single TCR/MHCp domain 
in the model. In contrast, the multifocal synapse of thymozytes is stable for 
hours. One reason for the pattern stability might be the thymozyte cytoskele- 
ton. Unlike the cytoskeleton of mature T cells, the cytoskeleton of thymozytes 
presumably remains in a mobile, nonpolarized state that still allows cell migra- 
tion [IT]. The few TCR/MHCp clusters of thymozytes may be coupled to the 
cytoskeleton, thus following its movements. 
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Appendices 

A. Continuum model for homogeneous membranes 

The membranes are regarded as thin elastic sheets that exhibit an average 
orientation parallel to a reference plane for which we choose Cartesian cood- 
inates x = {xi,X2). First, let us consider a single membrane. The shape of 
this membrane can be parametrized by the local displacement field h(x) which 
describes deformations from the planar reference state with h{x) = 0. The 
associated elastic energy has the general form 

noi{h}^nto{h} + nbc{h} (a.i) 

where the first term Titc represents the work against the membrane tension 
(Ti, which is conjugate to the total membrane area, and the second term Tibe 
corresponds to the bending energy which is governed by the bending rigidity 

Ki. 

The work against the tension ai is given by 

Htcih} ^ ai J d^x + (V/i)2 - 1 w y" i (V/i)2 (A.2) 

where the asymptotic equality holds to leading order in the gradients of h. This 
term has the same form as for a fluctuating interface. [H] The bending energy 
Tibe depends on the squared mean curvature (85j. When expressed in terms of 
the displacement field h, the mean curvature M is given by 

OAT- ^''^+(^20'(ai9l/^)-2(aiM(a2/l)(5l92/l) + (9l/l)2(5252/i) 

- [l + (V/i)2]3/2 ^^■■^> 

where di and 82 represent partial derivatives with respect to xi and X2, respec- 
tively. To leading order in the gradients of h, one has 

2M = [V^/i + 0{\'H^)]/[l + (V/i)2]3/2 ~ v^/i. (A.4) 

The bending energy is then given by 

Hbcih} = «i y" VTTiyW ~ j d'xl (V2/i)2. (A.5) 

where the asymptotic equality corresponds again to the small gradient limit. 
Thus, in this latter limit, the elastic energy Tic\ behaves as 

n^i{h)^ j (fx [i(Ti(V/i)V i/ti (V^/i)^] (A.6) 



^For a deformable membrane interacting with a planar surface, the separation field I as 
used in section [2. is identical with the displacement field h considered here. 
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which is identical with the expression ([2| if one sets h = I, ai = a, and ki = k. 

Next, let us consider two membranes which are, on average, both oriented 
parallel to the {xi,X2) plane. The two membranes are distinguished by the 
index j with j = 1,2. The membrane with index j is described by the displace- 
ment field hj = hj{x), which represents the local distance of membrane j from 
the reference plane, has bending rigidity Kj, and is subject to the tension aj. 
To leading order in the gradients of hi and /12, the elastic energy of the two 
membranes is then given by 

In addition, the two membranes interact via the effective potential Vmc which 
depends on the separation field 



l{x) = hi{x) - h2{x) > 



(A. 



where we used the convention that membrane 1 is located above membrane 2. 
The inequality l{x) > reflects the basic property that the two membranes 
cannot penetrate each other and, thus, exert a mutual hardwall potential. In 
general, the interaction potential V may depend both on the separation / = 
hi — /i2 and on the gradients of hi and /12. As before, we keep only the leading 
term in a gradient expansion and, thus, consider the simple interaction energy 



'H2An{hl, h2 



cPx 



V{hi - h2) = H2M{hl - h2}. 



(A.9) 



The effective Hamiltonian for the two membranes is then given by 

n2{huh2} - n2Mhi,h2} + H2M{hi ~ h2\ (A.IO) 



with 7i2,ei{/ii, ^2} as in (A.7 1 



We now change variables from the displacement fields hi and /12 to the 
separation field I ^ hi — h2 and another new displacement field m. The new 
field m is most conveniently described in Fourier space. Thus, let us define the 
wavenumber q = (91,92) and the Fourier transform / — f{q) of any function 

/ = f{x) via 

~f{q) ^ I d'x e-^^-^ix). (A.ll) 



The elastic Hamiltonian as given by (A.7 1 can now be expressed in terms of the 
Fourier transformed displacement fields hi — hi(q) and /12 = ^2(9) which leads 
to 

d'^q 

2 AiW;!"!! T 2 



'H2,c\{hi,h2} 



with 



(2^)2 



'2Xiiq)\hi\'+ '2X2(1)1^2 



a^q^ + Kjq'^. 



(A.12) 
(A.13) 
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The Fourier transform of the new displacement field I is simply given by 
I ~ hi ~ h2- The new displacement field m, on the other hand, is now defined 
via its Fourier transform 

m= ^1+ h2 (A.14) 

Xi + X2 Xi + X2 



with the g-dependent functions xi and X2 as given by (A. 13 1. 

When expressed in terms of the new displacement fields I and m, the effective 
Hamiltonian (A. 10) for two membranes is decomposed into two terms according 
to 

H2{hi,h2} = H{l} + 'H'{m} (A.15) 

and, thus, does not contain any crossterms involving both I and m. The effective 
Hamiltonian TL' {m\ does not depend on the interactions of the membranes and 
describes the diffusive motion of the displacement field m which is analogous to 
a 'center-of-mass' coordinate. 

The effective Hamiltonian 'H{V\, which governs the separation field I — hi ^ 
/i2, is given by 

WW = / (0 '^x{qW + Id'xV{l) (A.16) 
with the 'inverse propagator' 

X <? = -7 , X 2 I / , V^- A.17 

{(Ti + (72)q^ + (Kl + K2)q^ 

If the tensions and bending rigidities of the two membranes satisfy the linear 
relation 

ai/(T2 ^ ki/k2, (A. 18) 

the function x('?) simplifies and becomes 

2j1 + 2j2 Ki + K2 



After the inverse Fourier transform, the effective Hamiltonian (A.16) then has 
the simple form |45l |44j 



H{1} = J (Px [\ a {^If + i K (yHf + . (A.20) 
with the effective tension 

(7 = 0-10-2 /(c^i + cr2) (A.21) 

and the effective rigidity 

K = kiK2/(ki + K2) (A. 22) 

as in Q and ([5]). This form also applies (i) to tensionless membranes with 
fi = CT2 = 0, (ii) to very flexible membranes with ki = K2 = 0, (iii) to two 
identical membranes with gi = 0-2 and ki ~ K2, and (iv) to the limiting cases 
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in which one of the membranes becomes very tense and/or very rigid. In all of 
these cases, the effective Hamiltonian for two interacting membranes as given by 



(A. 20 1 is identical with the effective Hamiltonian for one dcformablc membrane 
interacting with another planar surface as described byQ - (|3| if the effective 
tension cricr2/(cri +(12) is identified with a and the effective rigidity kiK2/(ki + 
K2) with /t as in Q and ([5|, respectively. 

If all tensions and rigidities are finite but do not satisfy the relation as given 



by (A. 18 1, the function x{<l) has the asymptotic behavior 



xw ~ — - — q + -, — - — yT 1 (^-23 



for small q and 

for small tensions ai ~ (72- In both limits, the leading order term is again 



xiq) ~ — - — q + -7 — - — ^ q (A.24) 



contained in (A. 19) 



The effective tension a and the effective rigidity k define the length scale 
= (k/ct)^/^. On length scales that are large and compared to the elas- 
tic energy is dominated by the tension and rigidity term, respectively. In the 



tension-dominated regime, the effective Hamiltonian (A. 20 1 for two interacting 



membranes is identical with the effective Hamiltonian for wetting transitions 
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